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EXECUTIVE SUMMARY 


This Pinal Report (in five volumes) was prepared by 
Ocean Data Systems, Inc, for the Jet Propulsion Laboratory 
under Contract No. 954668 (as amended) in support of the 
SEASAT Program. 

The objective was to develop atmospheric analysis and 
prediction models of varying (grid) resolution, and to test 
the models using real observational data for the purpose of 
assessing the impact of grid resolution on short-range numerical 
weather prediction. The work statement was amended to include 
the performance of sensitivity tests using a coarse-mesh 
(63 K 63 X 5 level) prediction model in order to identify and 
order factors which might mask or impair the utility of SEASAT 
data on short-range weather prediction. Such factors included; 
initial conditions; topography; surface friction; latent heating; 
diffusion of moment\am and temperature; and computational devices 
such as tendency truncators, pressure smoothers, and temporal 
filters. 


A. Description 


ODSI has completed the sensitivity tests using model ver- 
sion PECHCV (63 X 63 x 5) . The model characteristics/parameters 
which were varied, together with forecast run descriptors, are 
listed below: 
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Temperature diffusion coef f i.::ient (K ) 

g 

= 10 (standard) ; Run SP-3 
= 10^; Run SP-4 

Momentum diffusion coefficient (K ) 

M' 

^14 = 10^ (standard) ; Run SP-3 

— Kjj = 10^; Run SP-2 

— = f (a) ; Run SP-1 

Surface drag coefficient (Cj^) 

Cq = 0.0015 (standard); Run SP-3 

— = 0.003; Run SP-7 

= 0.00015; Run SF-9 

Pres sure- smoothing coefficient (e) 
e =0.01 (standard); Run SP-3 
e = 0.02; Run SF-6 

— e =0,0; Run SP-5 

— e = f (t) ; Run SF-12 

Precipitation (large-scale plus convection) 

— On (standard) ; Run SP-3 

— Off; Run SF-14 


Temporal filtering 


— On (standard) ; Run SF-3 

— Off; Run SF-15 


Multiple modifications 

— Standard; Run SF-3 

— No diffusion, pressure smoothing, friction, radiation, 
condensation, evaporation, or sensible heating; 

Run SF-11 

— No radiation, evaporative/sensible exchanges at 
surface; Run SF-20 

— No diffusion (K^=Kj^=0 . 0) ; Run SF-10 


Initial wind specification 

— Analyzed (standard) ; Run SF-3 

— Derived from balance equation; Run SF-16 


Initial moisture specification 

— Parameterized using vorticity patterns (standard) ; 
Run SF-3 

— RH = 25%; Run SF-17 

— RH = 75%; Run SF-18 
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• Tendency truncation/Terrain tapering 

— Tendencies and terrain vanish south of 5®N, and 
increase smoothly with latitude to full value at 
20“N (standard) ; Run SP-3 

— Tendencies and terrain reduced by 50% south of 5°N, 
with latitudinal increases to full value by 20 °N; 

Run SF-19 

• Terrain specification 

— Tapered south of 20 "N to 5°N equatorward of which it 
vanishes (standard) ; Run SF-3 

Smooth earth; Run SF-13 

B. Findings 

These sensitivity studies have demonstrated that the 
effects of each process, procedure and parameter contributing 
to the computer forecast need to be understood, ordered and 
exploited. Many of the choices are model dependent. 

Physical effects are often small when compared to other 
sources of (possible) variation in a forecast. In fact, some 
computational/cosmetic devices can mask (or overwhelm) certain 
physical effects. Consider, for example, pressure-smoothing 
which has a greater impact on the forecast than surface friction. 
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mountains or vertical resolution. Momentum diffusion impacts 
on the forecast as much as (modeled) precipitation. The 
coefficients for both pressure smoothing and momentum diffusion 
should be specified only after careful evaluation of numerous 
forecasts. 

Table 1 summarizes some of the variations between two 
forecasts due to specified modifications in a procedure/ 
coefficient or effect. Precipitation amounts appear to be 
affected mainly by changes in the horizontal resolution and/or 
in the initial moisture specification. The average storm 
intensity is affected by horizontal resolution, too. Vertical 
resolution increase (from five to ten layers) had minimal 
impact on storm centers. Pressure smoothing and momentum 
diffusion have as much effect on lows as precipitation or 
surface friction. With respect to the kinetic energy change- 
predicted at the lowest model level, note that momentum 
diffusion has greater impact than tripling the horizontal 
resolution. Pressure smoothing has the same effect as does 
surface friction or precipitation. Both radiation and tem- 
poral filtering appeared to have only minor (if any) impact 
on a one-day forecast. 
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TABLE 1: VARIATIONS IN FORECASTS FROM SPECIFIED MODEL CHANGES. 


ITEM 

DIFFERENCE BETWEEN TWO FORECASTS 

PRECIPITATION 
AMOUNT (%) 

AVERAGE STORM 
INTENSITY (MBS ) 

KINETIC ENERGY 
CHANGE AT LOWEST 
LEVEL (%) 

RADIATION (ON-OFP) 

'V.l 

0 

0.0 

TEMPORAL FILTER (ON-OFF) 

1-2 

0 

» 

0.1 

FRICTION (2X CHANGE) 

2-3 

1.5 

14.6 

DIFFUSION (lOX CHANGE) 




- MOMENTUM 


3.0 

66.9 

- TEMPERATURE 

3-4 

1.0 

4.7 

PRESSURE SMOOTHING (2X CHANGE) 

4-5 

2.5 

14.6 

INITIAL WINDS (ANALY ZED-DERIVED) 

-v-S 

1.2 

12.7 

MOUNTAINS (IN-OUT) 

0,6 

1.1 

21.1 

PRECIPITATION (ON-OFF) 

N/A 

3.0 

14.0 

VERTICAL RESOLUTION (2X CHANGE) 

10-15 

0.5 

3.5 

HORIZONTAL RESOLUTION (3X CHANGE) 

80-100 

6.3 

57.7 

INITIAL MOISTURES (75%-25%) 

50-75 

2.6 

19.2 
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INTRODUCTION 


This Final Report (in five volxomes) was prepared by 
Ocean Data Systems, Inc. for the Jet Propulsion Laboratory 
under Contract No. 954668 (as amended) . 

The objective was to develop atmospheric analysis and 
prediction models of varying grid resolution, and to test the 
models using real data for the purpose of assessing the impact 
of resolution on short-range (1-2 days) forecasts. This work 
statement was amended to include the performance of model 
sensitivity studies, using a coarse-mesh (63 x 63 grid with 
five vertical layers) prediction model, in order to identify 
and to order some of the factors which might mask or impair 
the utility of SEASAT data. 

. This volume contains the descriptions and results of all 
sensitivity tests performed. Model parameters, procedures and 
effects were altered singly or in combination to determine the 
impact on short-range forecasts. These included such things as: 
initial conditions (winds; moistures) ; radiati.on; orography; 
surface friction; precipitation; diffusion (momentum; tempera- 
ture) ; tendency truncation; pressure smoothing; and temporal 
filtering. 


Numerical models, of course, are used to simulate the 
physical/dynamical processes, effects and interactions taking 
place in the real atmosphere. These models tend to differ 
from one another in terms of grid resolution, the physics, 
initialization, and numerical methods and procedures. Many 
of these procedures and parameters are discretionary in nature. 
Thus, there is a need to determine the impact of each procedure 
and parameter on the forecast. 


II. METHOD 


Since the inception of numerical weather prediction, 
numerous computational devices and procedures have been devised 
and implemented in order to improve either the efficiency, 
validity or stability of the solutions. Such devices have been 
modified or discarded with experience. The methods employed in 
NWP are such that smoothers, filters, tendency truncators, 
mountain tapering (and the like) are necessary to enhance the 
stability and appearance of the numerical forecasts. On the 
other hand, their repetitive use can (and does) have hainnful 
side effects. 

The purpose of this series of sensitivity tests with a 
coarse-mesh primitive-equation prediction model is to determine 
and to order the impacts of such procedures and devices , many 
of which are discretionary. The actual model (PECHCV) is 
described in detail in Volume III of this Final Report. The 
initial conditions for all tests are those for 12002, 22 April 
1976. Each of the tests will be introduced and described in 
subsequent paragraphs. 


A. Temperature Diffusion 


The diffusion terms of the thermodynamic energy equation 
is given by 

m^ . . 

D = D — ^ (7rV^T}y"i 

The diffusion coefficxent^ D/ varies with latitude 
according to the relation K (l-sin^<})) with 0.9 taken as a 
maximum for sin(J). K is taken to be 10®. 

purposes of this study / the model is executed 

with K = 10® (SP4). The notation (SF ) denotes the label 

by which this particular alteration will be referenced in 
the tables which follow. 

B. Momentum Diffusion 

The momentum diffusion terms of the momentum equations 
treated in the same manner as the temperature 
diffusion term above. For this study, runs were made where 
the coefficient of the momentum diffusion terms are: 10® 

(SF2) and varying values from 1.6x10® to 2.2x10® depending 
upon the level (SPl) . 
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c. 


Drag Coefficient 


Frictional dissipation is accounted for in the gross 
boundary layer (a = 1.‘0 - 0.8) by the stress terms Fx and 
Fy, respectively / of the momentum equations in the x and y 
directions. The friction in the x direction is given by: 


,n-l 

X 


(^) 

AaR '^s 't ’ 


n-1 


and in the y direction by 


,n-i 


^^n-l 

AaR ^ 


where = 0.3 [u^ + 


and Cq = 0.0015 for sea level and 0.0025 for non-sea level 
terrain in the baseline version. Cj^ is experimentally 
modified to .003 (SF7) and .00015 (SF9) for this study in two 
separate executions of the model. 


D. Pressure Smoother 

In order to control small-scale noise in the pressure 
field (PS ) f a non-linear smoothing operator is applied every 
time step. This smoothing operator has been described by 
Oliger and Wellck (1970) and is given by; 


APS 


dPs i 


max 




?Kh 


^^-■"»i-i,ii ("Si" - "®i-i,j> 


+ |pSi,j+i - ^§♦’■1 


- <") 


’i,j+l. 


where I 0Ps 1 is the maximum absolute difference between Ps 
at the point (i,j) and Ps at any of the four immediately 
neighboring points, and s is the smoothing coefficient which 
has been determined to be oi the order of 0.01 for control of 
computational noise. The smoothed surface pressure is then 
given by; 

Ps*^ (smoothed) = + APs 

where Ps^ denotes Ps at (i,j). For the purposes of this study, 
smoothing coefficients of 0 (SP5) , .02 (SP6) and values as a 
function of time (SF12) are demonstrated in three ^.uns. 
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E. 


Precipitation 


Both large-scale and convective-scale precipitation is 
modeled ordinarily. In the large-scale mechanism, 100% rela- 
tive humidity is the criterion. Layers beneath the precipi- 
tating layer are checked to see if evaporation is called for, 
before the precipitation reaches the ground. Three types of 
cumulus clouds are parameterized, two of which may precipitate. 
For our sensitivity runs, we turned off these mechanisms (SF-14) . 

F. Time Filter 

In order to achieve damping of the shorter wave length 
components to avoid non-linear instability and, in particular, 
to damp the computational mode, the simple leap-frog scheme 
is modified slightly to time filter the solutions of the 
momentum, thermodynamic and moisture equations. This procedure 
is known as time filtering and was originally used by Robert 
(1966) . 

Assume that the partial differential equation that we 
wish to approximate is 

H = 

Then, using the leap-frog difference scheme, we have 

= A^“^ + 2At F(A) 
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This step is modified by Robert to become 


= A^ ^ + 2 At F(A,A') 

a”^ = (1-a) (A')^ + j + (A')^"^^]/ 0<a<l 

Note that after A’ is computed at the n+1 time level, A' at 
the nth time level is converted to an A value by the mixing 
of time levels. This A is then used in the calculation at 
the n-1 time level. 

It can be shown that this procedure damps the meteoro- 
logical solution and heavily damps the computational mode 
that is present in the solution of the finite difference 
equations. For this study, the model is executed, in one 
case, without the assistance of the time filter (SF15) . 

G. Heating 

In the normal execution of the model, the contributions 
of solar radiation, terrestrial radiation and sensible heating 
and evaporation in the planetary boundary layer are computed 
once each hour and pro-rated over the next one-hour forecast 
period. For the purposes of this study, a selected run of 
the model is executed where these contributions are not 
operative (SF20) . 
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H. 


Initial Winds 


The initial winds in this model are normally obtained 
from an analysis procedure which gives analyzed winds on 
pressure surfaces from 1000 to 100 mb, and are computed 
geostrophically on the 50 mb pressure surface. The model 
uses wind data on eleven pressure surfaces; 950, 850, 700, 
500, 400, 300, 250, 200, 150, 100 and 50 mb. The initial 
winds on pressure surfaces are then vertically interpolated 
linearly in the logarithm of pressure to each model sigma 
surface. 

For the purposes of this study, the model is run once 
with the winds initialized by the solution of the divergence 
(balance) equation on constant pressure surfaces and the 
interpolation of the rotational component of the wind to 
sigma surfaces (SF16) . 

I, Moisture (Initial Relative Humidity) 

The model is affected by moisture, which is normally 
initialized by a parameterization method, in its lowest 
three levels. This initialization method is based on the 
intensity of the geostrophic relative vorticity which is 
used to define the initial moisture distributions on the 
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1,000, 900, 700, 500 and 300 mb pressure surfaces. These 
values are then interpolated linearly in the logarithm of 

pressure to the three sigma surfaces. 

in this study, two runs of the model were made in 
which a different constant value of EH was specified for the 
lower three levels. The value of 25% (SP17) was used for 
one run and 75% (SF18) for the other. 


Boundary Conditions 

The model is bounded in the horisontal by rigid, imper- 
meable vertical walls placed on the grid on the next-to- 
outermost row of grid points. The boundary is both insulated 
and slippery, that is, no heat or momentum exchange is 
permitted across the wall, but a parallel flow is permitted 

along the wall# 

There is also a provision in the model to truncate the 
computed tendencies south of a specified latitude. In the 
standard model (SF-3) , this procedure was operative. From 

to the pole, full values of the computed tendencies are 
used. Below 5°N, the tendencies are set to zero. In the 
region between 5“N and 20“N the amount of the computed 
tendencies used varies according to the relation [ (sin<j>-sin5 )/ 
(Gin20°-sin5“)l^. This procedure results in a persistence 
forecast in the equatorial areas, and acts as a computational 

buffer. 
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In this study, a special run (SF19) of the model was 
performed in which the truncation tendency was a linear value 
between 1 and .5 from 20°N to 5°N, and a constant .5 south 
of 5°N. The terrain is normally tapered south of 20°N by 
the tendency truncation factor, but in SF19 this tapering 
is omitted, 

K. Terrain 

The terrain height is normally area-averaged and gradient- 
limited. That is, the terrain gradient at any grid point is 
constrained to be less than 2000m per 381km, the basic 
63 X 63 grid size. In addition, the terrain is multiplied by 
the tendency truncator, thereby tapering it towards sea level 
around the edges of the grid. 

For this study, two runs of the model were made with 
specially initialized terrain. In one run, the terrain v/as 
set to zero {flat earth) (SF13) and in the other, the tapering 
procedure was omitted thereby representing terrain at its 
full smoothed height south of 20°N (see F19, above) . 

In the above expressions, 

u = X direction wind component 

V = y direction wind component 

TT = terrain pressure 
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T = temperature 

ra = map factor 

F^/Fy, = stress components 

= temperature diffusion term 
n = time step 


III. SYNOPTIC PATTERN 


All of the experiments described in the previous section 
were performed as twenty- four hour forecasts with initial 
conditions at 1200Z 22 April 1976 forecast to verify at 1200Z 
23 April 1976. The initial analysis at the surface is shown 
in Figure III-l and at 500 mb in Figure III-2. The verifying 
surface and 500 mb analyses are given in Figures I I 1-3 and 
III-4, respectively. 

It should be pointed out, once more, that these analyses 
are not the actual synoptic situations at the given times but 
a simulation of the atmosphere based upon the analysis tech- 
nique and its constraints, the available data and the quality 
of the "first guess" field. The verification analyses are in 
turn dependent upon the initial analyses since the forecasts 
from the initial analyses were constituents of the verifi- 
cation’s first guess fields. 

The 24-hour history of some of the primary systems of 
this synoptic pattern are given in Tables III-l and III-2. 

The changes in Figures III-5 and III-6 are due to a combin- 
ation of system development and movement. 

Tables have been prepared in this and the next section 
to numerically define the degree of error or change between 
two fields. The tables are sxunmarized by ten-degree latitude 
bands and for the entire northern hemisphere. The statistics 
that are summarized are: 


III-l 


Root mean square (RrciS) 


RI4S 


rp T 9 

E( F - •^V) 

N 


1/2 


Standard Deviation (SD) 


SD = 


- 

N 





1/2 


Mean = - "^V) 

N 


Tp = forecast parameter ' 

Ty = verifying parameter 
N = number of grid points 

as well as the algebraic maximum (I4AX) and minimum (IIIN) 
values. 

Table III-3 is the summary of the 24-hour changes for 
this synoptic period. The greatest change at the surface r 
associated with the East Asia low, was 40.8 mb. At 500 mb 
the greatest change was 381 m , associated with the same 
system. 
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FIGURE III-2: ANALYSIS, 500it±», 1200z 22 April 1976 

(Initial Condition) 
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FIGURE III-4: ANALYSIS, 500mb, 1200z 23 April 1976 

(Verification) 
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FIGURE III-6: ACTUAL CHANGE, 24-HOUR PERIOD, 500mb 

(60M Interval) 

Original page 
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TABLE III-l: SURFACE PATTERN HISTORY 1200Z 22 APRIL 

1976 to 1200Z 23 APRIL 1976. 



48N 

154W 

45N 

164E 

32N 

153E 

45N 

148E 

44N 

125E 

40N 

103E 


48N 

140W 

48N 

176W 

34N 

159E 

46N 

158E 

49N 

129E 

37N 

lOOE 

6 IN 

low 

47N 

37W 


47N 82W 


1004 


5 IN ^7W 
(split, 
low form- 
ing off 
Cape 

Hatteras) 
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TABLE III-2: 


500 MB PATTERN HISTORY 1200Z 22 APRIL 
1976 to 1200Z 23 APRIL 1976. 


LABEL 

TYPE 

INITIAL 
HT (M) 

INITIAL 

LOCATION 

FINAL HT (M) 
OR COMMENT 

FINAL 

LOCATION 

A 

Trough 

- 

along 118W 

slightly 

v;eaker 

along 114W 

B 

Trough 


along 170W 

low from 
Aleutians 

along 155W 

C 

Low 

5280 

48N 

117W 

5160 

48N 

126W 

D 

Low 

5340 

57N 

19E 

5340 

55N 

lOE 

E 

Low 

5280 

48N 

46W 

5280 

45N 

36W 

P 

1 

Trough 

- 

Wisconsin 
to North 
Carolina 

- 

i 

L. Erie to i 
SE US Coast j 



1200Z 22 APRIL 1976 TO 1200Z 
23 APRIL 1976 24-HOUR CHANGE 
STATISTICS . 


TABLE III-3: 


SURFACE 


FROM TO 
0 10 
n 20 
20 30 
30 40 

40 50 

50 60 

60 70 

70 50 

80 90 

0 90 


FROM TO 
0 10 
10 20 
20 30 

30 40 

40 50 

50 60 

60 70 

70 80 

80 93 

0 90 


RMS 
.141E+01 
.168E+01 
.260E+01 
.515E+01 
♦107E+02 
,10SE+02 
.743E+01 
, 883E+01 
.614E+01 
.514E+01 


RMS 

.157E+02 

.211E+02 

.29EE+02 

.661E+02 

,101E+03 

.912E+02 

.715E+02 

.717E+02 

.6Z5E+02 

.503E+32 


SO 

.117E+01 
. 1686+01 
-,25SE+01 
.515E+01 
.105E+02 
.108E+02 
,743E+01 
.859E+01 
, 373E+01 
, 514E+01 


500 

SD 

.119E+02 
. 206E+02 
,297E+02 
.655E+02 
.lOlE+03 
.910E+02 
.715E+02 
.713E+02 
.618E+02 
, 503E+02 


MEAN 
,791E+00 
-.125E-01 
.345E+00 
, 409E-02 
-,200E+01 
-*375E+00 
. 726E-01 
.205E+01 
,487E+01 
, 179E+00 


MB 

MEAN 
, 103E+02 
. 462E+01 
-.112E+01 
-♦861E+01 
-» 546E + 01 
. 559E+01 
, 212E+01 
-.765E+01 
894E+01 
. 260E+01 


MAX 

•499E+01 

.445E+01 

.963E+01 

.161E+02 

.210E+02 

.210E+02 

,128E+02 

.144E+02 

.120E+02 

.210E+02 


MAX 
. 439E + 02 
,550E+O2 
.113E+03 
.150E+03 
.2Q7E+03 
,2186+03 
.193E+03 
.120E+Q3 
,114E+03 
.2186+03 


OgfGlNAi; PAQ& Xi: 


MIN 
.242E+01 
» 4 5 5 E +0 1 
, 835E+01 
. 195E+02 
.408E+02 
,26BE+02 
,2406+02 
.2346+02 
♦623E+01 
.4086+02 


HIN 
.292E+02 
.423E+C2 
.1526+03 
. 3156+03 
. 381E+03 
.2586+03 
.1886+03 
.1456+03 
. 1096+03 
.3816+03 


IV. EVALUATION OP FORECASTS 
A. General 

The run identified as SF3 (see Figures Iv-1 and IV-2) 
is the baseline version of the model for this study. That 
xs, all of the engineering devices and coefficients were 
set to their generally accepted operational values. The 
other runs xn the study are compared to this version as 
well as the verification analysis (error charts) and in 
some cases to each other. 

Analyzing the performance of SP3 in detail and with 
the aid of its error charts. Figures IV-3 and IV-4, it is 

seen that the following characteristics were demonstrated by 
the model; 

Surface - Low A moved too slowly and deepened in place 
for an error of 12 mb. Low B showed good movement but was a 
bxt too slow for an error of 8 mb. Most of the development 
and movement of low C was picked up but it tended to erode 
the mid-Pacific high too rapidly for an error of 10 mb at 
about 170E. 

The model performed very well in predicting the deepen- 
ing and movement of low E, Due to the tight gradient of the 
low and slightly too much eastward extension, a relatively 


mxnor error of about 10 mbs Xi7as produced. This compares 
favorably with the actual change of 40 mb in this area. 

The forecast position, intensity and central pressure 
of low H were excellent and resulted in almost no error. 

The forecast status of low I was also excellent and 
included the low development east of Cape Hatteras. 

500 MB - Figure IV— 4 shows little error of greater 
than SOM except for the following: 

In the Gulf of Alaska, the model developed the strong 
westerly flow into a slight ridge rather than a weak trough 
in advance of the approaching low which resulted in the 
largest error of 150+M. This error was very close to the 
actual 24-hour change of 18 OM. 

The intensity and deepening of the western Asian low 
were excellent. There was, however, a bit too much eastward 
extension of the low circulation (reflected at the surface) 
which lead to a 90M error. This was an excellent forecast 
in an area of a 381M actual change. 

The model attempted to develop a weak low north of 
Thailand for an error of 12 OM. 

The error summary for SP3, Table IV-1, shows the largest 
error at the surface to be 12.7 mb and 163M at 500 mb. 
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TABLE IV-1: ERROR STATISTICAL SUMMARY, 

SP3 (BASELINE RUN) . 


FRDH TO 
0 10 
10 20 
20 30 

30 40 

40 50 

50 60 

60 70 
70 00 

80 90 

0 90 


FROM TO 
1 10 
10 20 
20 30 

30 40 

40 50 

50 60 

60 70 

70 £0 
80 90 

0 90 


RHS 

.710E+00 

.116e+01 

♦292E+01 
.339E+01 
.380E+01 
. 434E+01 
,332E+01 
.373E+01 
.239E+01 
.250E+01 


RHS 
. 165E+02 
.266E+D2 
.228E402 
.366E+02 
.476E+02 
.3468+02 
.251E+02 
.402E+02 
.439E+02 
,2S5E+02 


SURFACE 

SO 

.672E+00 
.116E+01 
.279E+01 
.304E+01 
.321E+01 
.346E+01 
. 275E + 01 
.357E+01 
.151E+01 
. 236E+01 


500 MB 
SO 

.147E+02 

.247E+02 

.190E+02 

.340E+92 

.476E+02 

,3476+02 

.236E+02 

.226E+02 

,117E+02 

,2776+02 


MEAN 
,23l)E+00 
, 185E-01 
. 860E+00 
.150E+01 
.204E+C1 
.260E+01 
.185E+01 
.1C9E+01 
.164E+01 
,819E+00 


MEAN 
.757E+C1 
,991E+01 
. 127E+02 
.137E+02 
, 420E+O1 
“,2e3E+01 
-.876E+01 
331E+02 
--.423E + 02 
•665E+01 


MAX 

.242E+01 

.412E+01 

.926E+01 

,112E+02 

.127E+02 

.122E+02 

,7756+01 

,894E+C1 

,5736+01 

,127E+02 


MAX 
.5446+02 
.6236+02 
.6376+02 
.1276+03 
. 1126+03 
,7776+02 
,5386+02 
.9116+01 
-.2346+02 
.1276+03 


HIN 
155E+01 
-.2826+01 
-.5626+01 
-.5676+01 
-.7746+01 
-. 1016+92 
-.7256+01 
-.9006+01 
-.7736+00 
-.1016+02 


HIN 

-.3526+02 

-.4966+02 

-.4236+92 

-.1146+03 

-.1636+03 

-.1216+03 

-.8656+92 

-.1046+03 

-.6586+02 

-.1636+03 
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B. Variable Momentum Diffusion (SFl) 

In this run, the momentum diffusion coefficient was 
varied as a function of level as follows: 

Level Coefficient 

1 1.0x10® 

2 1.9x10® 

3 2.2x10® 

4 2.2x10® 

5 1.0x10® 

Figure IV-5 shows that the SFl run closely resembles 
the baseline at the surface with generally minor differences 
except for the eastern Asia low. There was better control of 
the eastward extension of this low, and consequently the 
error chart, Figure IV-6, shows a far less complex error 
pattern in this area. Figure IV-7, the difference between 
SFl and SF3 at the surface, shows as much as a 10-mb difference 
here. The statistical summary of SFl, Table IV-2, has an 
overall lower maximum error at the surface than seen in 
the baseline, Table IV-1. 

The differences at 500 mb are less pronounced than at 
the surface. Figure IV-8 is the SFl 500 mb forecast. Figure 
IV~9 is the SFl forecast error, and Figure IV-10 is the 
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difference between SPl and SP3 at 500 mb. The latter two 
figures show an improved treatment of the east Asia low by 
SPl. Table IV-2 shows very little departure from the base- 
line in the overall error statistics at 500 mb. 


















FIGURE IV 7: DIFFERENCE BETWEEN BASELINE AND SFl, SURFACE 
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TABLE IV-2: ERROR STATISTICAL SUMMARY, SFl 

(MOMENTUM DIFFUSION COEFFICIENT 
FUNCTION OP TIME) . 


SURFACE 


FROM TO 
0 10 
10 20 
20 3 0 

30 <!t0 

AO 50 
50 60 

60 70 
70 QO 
60 90 

0 90 


FROM TO 
0 10 
‘ 10 20 

20 30 

'30 AO' 
AO 50 
- 50 60 

60 70 
70 80 
00 90 
0 90 


RMS 
,711E+00 
. ,116E+01 
.293E+01 
.327E+01 
. 353E+01 
.395E+01 
.32AE+01 
.372E+01 
, 233r:+0a 
,2AOE+01 


RMS 
,166E+02 
.269E+02 
.217E + 02. 
.355E+02 
.50AE+02 
, , 357E + 02 
.2A5E+02 
.A07E+02. 
. AA6E + 02 
.238E+02 


SO 

. 669E+00 
,116E+01 
.270E+01 
.298E+01 
.335E+D1 
. 321E+01 
.270E+01 
.3586+01 
.159E+01 
.228E+01 


500 

SD 

.lA7£+02 
.2A7E+0? 
.175E+02 
.332E+02 
.50AE+02 
. 355E+02 
.230E+02 
.233E+02 
. 112E+02 
.281E+02 


MEAN 
.2ALE+3Q 
.S05E-01 
. 8A7F + CO 
.13AE+01 
.1A2E+01 
.229E+C1 
. 179E+01 
.1*03E + 01 
•171E+01 
.7A0E+0O 


MB 

MEAN 
.766E+01 
.106E+C2 
.126E+02 
. 125E+02 
238E+00 
“.367E + 01. 
8A8E+01 
33AE+02 
-.A3AE+02 
.627E+01 


MAX 
. 2A2E+01 
.A22E+C1 
. 969E+01 
.117E+02 
.953E+01 
.963E+C1 
.761E+01 
.859E+01 
.53AE+01 
, 117E+02 


MAX 
.539E+02 
•630E+02 
,638£+D2 
.llAE+03 
.113E+03 
.625E+02 
.A92E+02 
.725E+01 
— • 2 6 6 E + 02 
, llAE+03 


MIN 

”, 156E+01 
-.276E+01 
-.615E+01 
58AE+01 
-.8A7E+0X 
77AE+01 
7A6E+01 
-.933E+01 
-.1A2E+0I 
-.933E+01 


MIN 

352E+02 

-.A97E+C2 

-.339E+02 

-.113E+03 

165E+03 

-.1236+03 

-.6A6E+02, 

-. 1066 + 03 : 

-.68AE+02 

165E+03 


C. Momentum Diffusion Coefficient = 10^ (SF2) 

By reducing the momentum diffusion coefficient by one 
order of magnitude, the importance of this device is seen in 
the SF2 surface and 500 mb forecasts. Figures IV-11 and IV-12, 
respectively. Most of the synoptic systems are in about the 
same locations and have about the same intensities as the 
baseline case, but considerable instability "noise" has 
developed. It is likely that this version would become compu- 
tationally unstable before 72 hours of forecast time. 

Systems that appeared to be developing only a little 
too rapidly on the baseline run are developing much more 
rapidly on SF2. For example, on the baseline run surface 
chart, there is a hint of a false extension of the east Asia 
low to the northeast that lead to an error of about 10 mb. 

SF2 develops this into a deep, intensive low west of Kamchatka 
and the main body of the low is split into two distinct 
centers with southern extensions. At 500 mb, in addition to 
corresponding incorrect development, the neophyte low north 
of Cambodia on the baseline has become an intense little low 
on the SF2 . 

The corresponding error charts. Figures IV-13 and IV-14, 
show patterns quite similar to the baseline, except for greatly 
increased errors in the regions discussed above. In fact. 
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Figures IV-15 and IV-16 show that the error differences are 
almost entirely confined to these areas. The statistical 
summary of SP2, Table IV-3, shows,, among other facts, that 
the size of the maximum errors are about double those of the 
baseline. 

The problem pointed out here as well as the slightly 
reduced errors in the same areas seen on the SFl run might 
suggest more favorable forecasts with a slight increase in 

g 

the momentum diffusion coefficient to possibly 1.5x10 or 
2.0x10^ vice 1.0x10^. 

One should take note of the radical changes in the 
forecast which can be produced by altering the momentum 
diffusion coefficient. 
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IV-12: FORECAST, SOOitita, SF2 
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FIGURE IV-14: ERROR PATTERN, 500mb, SP2 
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FIGURE IV-15 ; DIFFERENCE BETWEEN BASELINE AND SF2 , SURFACE 
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;j FIGURE IV-16 1 DIFFERENCE BETWEEN BASELINE AND SF2^ SOOmb 
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TABLE IV“3: ERROR STATISTICAL SUMMARY, SP2 

(MOMENTUM DIFFUSION COEFFICIENT 
= 105 ) . 


SURFACE 


FROM TO 
0 10 
10 20 
20 30 

30 40 

40 50 

50 60 

50 70 

70 80 

• 60 90. 

0 90 


FROM TO 
0 10 
10 20 
20 30 

3J 40 
40 50 

50 60 
50 70 

70 80 

80 90 

0 90 


RMS 
.708E+00 
.122E+01 
.326E+01 
. 426'c+Ol 
.498E+01 
«566E+01 
.357E+01 
.384E+01 
e254E+01 
.302E+01 


PHS 
♦164E+02 
.268E+02 
.276E+02 
.444E+02 
.522E+02 
.410E+02 
.271E+02 
.406E+02 
, 425E+02 
.315E+02 


SD 

.676E+00 
,122E+01 
.309E+01 
.383E+01 
.410E+01 
.490E+01 
. 301E+01 
,36BE+0l 
. 149E+01 
.287E+01 


500 
S 0 

.145E+02 
.252E+02 
.234E+02 
.414E+02 
.515b+02 
. 409E+02 
.257E+D2 
,232E+D2 
.131E+02 
, 306E+02 


MEAN 
.212E+00 
~.402E“01 
.106E+Q1 
. 186E+01 
•283E+01 
,284E+01 
. 193E-5-01 
.llOE+01 
,206E+01 
.959E+0D 


MB 

MEAN 
.741E+01 
.912E4-01 
. 147E+02 
. 162E+02 
•881E+01 
-.260E+01 
~.898E+01 
-.334E+02 
-.404E+02 
. 743E+01 


MAX 

•244E+01 

.394E+01 

.104E+02 

.158E+02 

,162E+02 

.265E+02 

.B49E+01 

.919E+01 

.613E+01 

.265E+02 


MAX 

.5575+02 

.6315+02 

.875E+02 

.214E+03 

.136E+C3 

.2075+03 

.592E+02 

.102E+02 

-.191E+02 

.214E+03 


HIM 

“.154E+01 
-".291E + 01 
-.567E+01 
~,958E+01 
-.993E+01 
-.125t+02 
-.747E+01 
“,659E+ni 
-.229E+00 
“.125E+02 


HIM 
352E+02 
-.490E+02 
".553E+02 
-. 110E+C3 
“. 157E+03 
-.117E+03 
-.898E+D2 
-.105E+03 
“» 666E + 02 
157E+03 
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D. Temperature Diffusion Coefficient = 10 (SF4) 

An examination of the SP4 surface and 500 mb levels/ 
Figures IV-17 and IV-18/ shows only a hint of the noise 
seen in SF2. The error charts/ Figures IV-19 and IV'-20/ are 
only modestly different from those for the baseline which is 
clearly indicated by the lack of major difference centers on 
Figures IV-21 and IV-22. Table IV-4 shows that the statistical 
summary for SF4 is not significantly different from that of 
the baseline. 

In general, a reduction of the temperature diffusion 
coefficient by one order of magnitude did not have an impor— 
tant effect upon this 24— hour forecast. That is not to say 
that a longer run, a different synoptic situation or a 
further reduction of the coefficient would not show a 
deleterious effect upon the forecast. 
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FIGURE IV-18: FORECAST^ 500inbr SF4 
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FIGURE IV- 22: DIFFERENCE BETWEEN BASELINE AND SF4, 500mb 
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TABLE IV-4: ERROR STATISTICAL SUMMARY, SF4 

(TEMPERATURE DIFFUSION COEFFICIENT 
= 105 ) . 


SURFACE 


FROM -‘TO . RHS SD 

® 10 ■»713E+0D .675E+00 

10 .20 '.115E+0X .115E+01 

20 30 , 4308E+01 *289E+01 

30 40.-,365E+01 .336E+01 

40 50 .383E+01 ,319E + r)l 

50 60 .439E+01 »347E+Q1 

6-0 70 .338E+01 .274B + 01 

70 .80 .337E+01 .324E+01 

80 90 .177E+01 .139E+01 

9 ^90 .257E+01 .242E+01 

♦ ' * * ■ 


MEAN 

MAX 

MIN] 

.230E+00 

.244E+01 

-.156E+01 ' 

. 804E-01 

.407E+01 

-.284E+01 : 

• 108 E+0 1 

.980E+01 

-.615E+01 1 

• 143E+01 

.124E+02 

716E + 01 

• 212E+01 

.142E+02 

3Q8E+01 

.266E+01 

.146E+02 

890E+Q1 

.197E+01 

*880E+01 

733E+01 

.930E+00 

,643E+01 

818E+01 

•llOE+Ol 

.394E+01 

~.193E+01 

.867E+00 

• f 146E + 02 

890E+01 


500 MB 


from to phs so 

0 10 .165E+0H .147E+02 

10 .20 .266E+02 .250E+02 

20 30 .223E+02 .193E+Q2 

30 . 40 .356E+02 .336E-h02 

40 -50 .476E+02 .471E+02 

5Q 60 ,342E+C2 *342E+02 

60 70 .253E+02 .241E+02 

70 80 .AOCE+02 .243E+02 

80 -90 .377E+02 .146E+02 

0 90 .282E+02 .274E+02 


MAX MIN 

#750E+G1 4543E+02 -.351E+02 

.923E+01 .618E+02 -.^94E+02 

#113E+02 »619E+02 -.437E+02 

«118t+02 .142Ef03 -.119E+03 

.689E+01 .112E+03 -.166E+03 

-.176E+00 .795E+02 -.123E+03 

-.767E+01 .604E+02 “,910E+02 

-.318E+02 ■165E+02 -.106E+03 

-.3^i8E + 02 -.107E+02 -,611E+02 
.657E+01 .142E+03 -.166E+03 
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E. Smoothing Coefficient = 0 (SF5) 

The smoothing operation is performed on the surface 
pressure;, therefore, the effect of eliminating it produced 
less than SOM of change at 500 mb in SF5. 

The surface analysis. Figure IV-23, is almost identical 
to the baseline, and as seen in Figure IV-24, there is little 
more than 4 mb difference betv/een the baseline and SP5. 

The error summary {Table IV-5) is quite similar to that 
of the baseline. 

The results of omitting the smoothing function are 
cumulative and the errors that would be introduced accelerate 
and spread upward through the levels with time so that by 
72 hours there would be severe instability throughout the 


forecast domain. 










TABLE IV~5: ERROR STATISTICAL SUMMARY;. SF5 

(SMOOTHING COEFFICIENT = 0) . 


SURFACE 


FROM 

TO 

RMS 

SD 

MEAN 

MAX 

MIN 

0 

10 

.714E+00 

.680E-i-00 

.2I5E+00 

.247E+01 

-♦ISTE+Ol 

10 

20 

.119E+01 

.119E+01 

1^3E-01 

.406Ei01 

-.2S9E+01 

20 

30 

,293E+01 

.26^E+01 

.706E+00 

.943E+01 

--.5906 + 01 

30 

^0 

.344E+01 

.318E+01 

.131E+01 

.111E+C2 

“.727E+01 

^0 

50 

.393E+01 

.331E+01 

.220E+01 

.149E+02 

-,822E+01 

50 

60 

.A55E+01 

.377E+01 

.255E+01 

.1^(4E + C2 

-.1106+02 

60 

: 70 

.321E+01 

.282E+01 

. 152E+01 

.ai5E+01 

-.7766+01 

70 

80 

. 37^E+01 

.364E+01 

.829E+00 

.8e6E+01 

-.8536+01 

80 

^0 

,153E+01 

.127L+01 

.B51E+00 

.380E+01 

-.117E+01 

0 

90 

.255E+01 

.244E+01 

.746E+00 

,149E+02 

-.1106+02 





500 

MB 



FROM 

TO 

RMS ' 

SD 

MEAN 

MAX 

MIN I 

0 

10 

.1636+02 

, 1476+02 

, 7Z1E+01 

.5416+02 

-.354E+02 ' 

10 

20 

.267E+02 

.248E+02 

.9975+01 

,6286+02 

-.4966+02 

20 

30 

,2316+02 

.1906+02 

,3306+02 

,6326+02 

-.416E+02 

30 

; 40 

.380E+02 

,3516+02 

,1486+02 

.133E+03 

-.11,46+03 i 

40 

5D 

, 435E + 02 

.482E+02 

,5196+01 

.1156+03 

-,162E+03l 

50 

60 

, 358E+02 

.357E+02 

-.2316+01 

.8416+02 

-.123E+03 

60 

70 

.2536+02 

.2366+02 

-.891E+C1 

.5296+02 

-,8686+02 

70 

00 

.417E+02 

.237E+G2 

-.3436+02 

.7B7E+01 

-.1076+03 

60 

90 

.455E+02 

.1226+02 

-.4396+02 

-.2596+02 

-.6806+02 

0 

90 

.290E+02 

.282E+02 

, 680E+01 

,1336+03 

-.1626+03 


;l; 

T 

I 
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F. 


Smoothing Coefficient = .02 (SF6) 


As in SF5, the effect of altering the smoothing co- 
efficient is generally undetected at 500 mbs. At the surface. 
Figure IV-25, there is little obvious difference between SP5 
and SP6 except that in the latter case the lows are generally 
not as deep and the gradients around the dominant lows are 
weaker. The intensity of the differences of SF6 from the 
baseline, in Figure IV-26, are generally of the same magnitude 
as seen in Figure IV-24, but of the opposite sign. The total 
difference between SF5 and SF6 therefore is as much as 10 mb, 
shown in Figure IV-27. 

Table IV-6, the statistical summary for SP6, has little 
overall difference from the baseline stimmary but a slight 
shift of absolute maximum values to the negative. 
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FIGUEE IV-27; DIFFERENCE BETWEEN SF5 AND SF6 , SURFACE 
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TABLE IV- 6: ERROR STATISTIC2U^ SUMMARY, SP6 

(SMOOTHING COEFFICIENT = .02). 


SURFACE 


SD 

• 696 E +00 
.115E+01 
.278E+01 
.296E+01 
.3^^E+01 
.359E+01 

, 2 a 9 E +01 

.398E+01 

.212E+01 

.2't2E+01 


MEAN 
,181E+00 
-,8A5E-02 
. 929E+00 
, 158E+01 
. 179E+01 
,258E+01 
. 207E+01 
.127E+01 
, 260E+01 
.816E+00 


500 MB 


FROM 

TfJ 

0 

10 

10 

20 

20 

30 

30 

40 

40 

50 

50 

60 

60 

70 

70 

80 

80 • 

90 

0 

90 


FROM TO 
0 10 
10 20 
20 30 

30 AO 
AO 50 
50 60 

60 70 

70 60 

80 90 

0 90 


RMS 
.719E+00 
. U5 n+01 
.293E+01 
.335E+01 
.388E+01 
* AA2E+01 
. 356E+01 
.A18E+01 
.335E+01 
.255E+01 


RMS 

,166E+02 

.266E+02 

,225£+02 

.35AE+02 

.A73E+Q2 

„3A0C-+02 

.252E+02 

,392E+02 

*A26E+D2 

.281E+02 


SD 

. 1A7E+02 
.2A7E+02 
, 190E+02 
.331E+02 
.A71E+02 
, 3366+02 
s236E+02 
.223F+02 
.1156+02 
.27AE+02 


MEAN 
.757E+01 
.973E+01 
.1216+02 
.1256+02 
.3A1E+01 
327E+01 
-.869E+01 
-.3226+02 
-. AlOE+02 
.6326+01 


MAX 
.236E+01 
.AUE + 01 
.91AE+01 
.1106+02 
.106E+02 
.116E+02 
.911E+01 
,1086+02 
.7386+01 
.116E+02 


MAX 
o5A5E+02 
.62AE+02 
.6A1E+02 
,1226+03 
.nOE + 03 
.713E+02 
. 5A 6E+02 
.960E+01 
-.2166+02 
.122E+03 


HIN 

160E+01 

-.28AE+01 

-.6A8E+01 

-,550E+ni 

-.103E+02 

102E+02 

-.7236+01 

-.9926+01 

-.1106+01 

-.1036+02 


HIN 

-. 3536 + 02 
4985+02 
-.4286+02 
-.1156+03 
-.1646+03 
-.1196+03 
-.8586+02 
-.1016+03 
-.6356+02 
-.1646+03 
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G. Smoothing Coefficients as a Function of Time (SF12) 

In SP12, the smoothing coefficient was increased each 
time step linearly from 0 initially to .01 at 24 hours. The 
output products of this run showed no perceptible difference 
from the baseline run, and the statistical simimary (Table 
IV-7) showed only slight variations from the statistics pro- 
duced by the baseline. 
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TABLE IV- 7: 


ERROR STATISTICAL SUMMARY, SP12 
(SMOOTHING COEFFICIENT = FUNCTION 
OF TIME) . 


SURFACE 


FROM TO 
0 10 
10 20 
20 BO 
30 ^0 

40 50 

50 60 

60 70 

70 80 

80 00 
0 90 


RMS 
.712E+00 
. 118E+01 
.292E+01 
.342E+01 
.387E-i-01 
.441E+01 
.325E+01 
.354E+01 
.194E+01 
.252E+01 


SD 

. 572E+00 
.118E+ni 
. 281E-f-01 
. 310E+01 
. 321E+01 
.354E+01 
. 274E+01 
.351E+01 
. 134E+31 
.238E+01 


MEAN 
. 236E+00 
. 253E-01 
.817E+00 
. 146E+01 
. 217E+01 
. 264E+01 
. 174E+01 
» 998E + 00 
.140E+01 
. 815E+00 


MAX 
.245E+01 
.415E+01 
.939E-f01 
..112E+02 
.136E+02 
.130E+02 
. 767E+01 
.845E+01 
.493E+01 
.136E+D2 


NINi 
154E+01 i 
-.203EfOl 
579E+01 i 
-.638E+01 . 
781E+01 I 
102E+02 
740E+01 
-.856E + 01 i 
-.594E + 00 ! 
102E+Q2 i 


FROM TO 
0 10 
10 20 
20 30 

30 ^0 
40 50 

51? 60 

60 70 

70 80 

00 90 

) 90 


RMS 
. 165E+02 
.266E+02 
.229E+02 
*373E+02 
. 482E+02 
*353F+Q2 
.252E+02 
.411Ef 02 
.448E+02 
. 208E+O2 


500 

SO 

. 147E+02 
. 247E4-02 
. 190E+92 
. 345E+0Z 
.479E+02 
, 352E+32 
. 236E+0H 
, 234F+02 
. 120E+02 
.260E+02 


MB 

MEAN 
»743E+01 
. 994E+01 
.129E+Q2 
. 143E+02 
,483E+C1 
257E+01 
891E+01 
“•* 338E+02 
-.431E+C2 
. 675E+01 


MAX 
. 543E+02 
.623E+02 
.636E+02 
.130E+03 
.114E+03 
.779E+02 
.532E+02 
.076E+01 
-.246E+02 
.130E+03 


MIN 
353E+02 
495E+02 
-.419E+02 
114E+03 
163E+03 
-.122E+13 
-.871E+02 I 
-.106E+03 1 
-.675E+02 ' 
163E+03 
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H. Time Filter Removed (SP15) 

The removal of the time filter produces slight phase 
errors that increase non-linear ly with time. In the first 
24 hours, however, there was little to observe in the way 
of inconsistencies with the baseline. The statistical 
summary of errors for SF15, Table IV-8, confirms that at 
this point, SF15 is very much like the baseline. 
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ERROR STATISTICAL SUI4MARY, SP15 
(NO TIME FILTER) . 


TABLE IV-8j 


SURFACE 


FROM TO 
0 10 
10 20 
20 30 

30 

^0 50 

50 60 

60 70 
70 90 

80 90 

0 90 


FROM TO 
0 10 
1C 20 
20 30 

30 40 

^0 50 

50 60 

60 70 

70 PO 
80 90 

0 90 


RMS 
.710E+00 
,117E+0i 
.292E+01 
, 341E+01 
.383E+01 
. 438E+01 
.333E+01 
.371E+01 
,2396+01 
.251E+01 


RMS 
•165E+02 
..266E + 02 
.227E+02 
.366E+02 
.4786+02 
.349E+02 
.254E+02 
.406E+02 
,44^E+02 
.285E+02 


SO 

.672E+0P 
, 117E+01 
.278E+01 
.304E+01 
,322E+Q1 
.351E+01 
, 276E+01 
.355E+01 
.154E+01 
.237E+01 


MEAN 
. Z31E+00 
.416E-01 
, 89QE+00 
.154E+01 
.207E+01 
,263 E+01 
,186E+01 
, 109E+01 
.163E+01 
, 839E+00 


MAX 
.243E+01 
.416E+01 
.932E+01 
,112E+02 
.127E+D2 
, 123E+02 
,76BE+01 
.680E+01 
.583E+01 
.127E+02 


MAX 

,5^4E+02 

.622E+02 

.627E+02 

.124E+03 

.114E+03 

.V87E+02 

.544E+02 

.939E+01 

-.243E+02 

.124E+03 


MIN 

-.154E+01 
-,278E+01 
-.588E+01 
-♦587E+01 
-.777E+01 
lOlE+02 
716E+01 . 
-,896E+0li 
-,799E+00i 
lOlE+02! 


MIN 

-.353E+02 
-,^97E+02 
-.421E+02 
-♦115E+03 
-.164E+03 
120E+03 
~.E81E+02 
-,i04E+03 
676E + 02 
-.164E+03 


500 MB 


SD 

.147E+02 
.247E+02 
.190E+02 
.340E+02 
.477E+02 
.347E+02 
.237E+02 
. 228E+02 
. 120E+02 
,278E+02 


MEAN 
.754E+01 
.978E+01 
.124E+02 
. 136E+02 
.395E+01 
318E+01 
-.907E+01 
-.335E+02 
“.^27E+02 
. 648E+01 
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I. Drag Coefficient = .003 (SF7) 

The frictional drag device acts primarily on the lowest 
level and eventually affects the upper levels. By doubling 
the drag coefficient to .003, little effect is seen at the 
surface or 500 mb in 24 hours. The departure from the base- 
line at the surface. Figure IV-28, shows a maximum of about 
-4 mb north of Korea. Figure IV-29, the 500 mb level, shows a 
departure of only 30M. The statistical error summary for 
SP7, Table IV-9, closely resembles that of the baseline with 
a slight shift of the maximum value to the negative. 
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FIGURE IV-28: DIFFERENCE BETWEEN BASELINE AND SF7, SURFACE 
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FIGURE IV-29: DIFFERENCE BETWEEN BASELINE AND SF7 , SOOmb 





TABLE IV- 9: 


ERROR STATISTICAL SUMMARY, SP7 
(DRAG COEFFICIENT = .003). 


SURFACE 


FROM 

TO 

0 • 

10 

10 

20 

20 

30 

30 

40 

40 

50 

50 

60 

60 

70 

70 

80 

80 

90 

0 

90 


RHS 
.710F+00 
.116E+01 
.298E401 
.333E+01 
. 378E+01 
*'f27E+01 
•33 lE+01 
•386E+01 
•241E+01 
•250E+01 


SO 

.670E+Q0 
.116E+01 
.280E+01 
.2996+01 
.348E+01 
.364E+ni 
. 280E+01 
.374E+01 
* 165 E+0 1 
.237E+01 


MEAN 
•237E+00 
•742E-31 
. 102E+01 
•145F+01 
.143E+01 
.222E+01 
. 175E+01 
. 954E+0D 
. 176E+01 
.775E+00 


MAX 
•242E+01 
.4256+01 
.9206+01 
.113E+02 
. 115E+02 
•llOE+02 
•763E+01 
.919E+01 
•583E+01 
. 116E+02 


FROM 

TD 

0 

10 

10 

20 

, 20 

30 

30 

40 

40 

50 

50 

50 

60 

70 

70 

80 

80 

90 

0 

90 


RMS 

.1666+02 
. 267E+02 
.2296+02 
.3596+02 
.4686+02 
. 363E+02 
.255E+02 
.419E+02 
.4476+02 
.288E+02 


500 

SO 

.1476+02 
.2476+02 
. 190E+02 
.3386+02 
. 488E+02 
.3596+02 
.236E+02 
.2406+02 
.123E+02 
.282E+02 


MB 

MEAN 
. 763E+C1 
.lOlE+02 
. 129E+02 
0 122E+02 
.1706-01 
. 541E+01 
“ . 96 8 E + 0 1 
-. 343E+02 
-.430E+02 
. 599E+D1 


MAX 
.546E+02 
.625E+02 
.6416+02 
.1256+03 
.1106+03 
.718E+02 
.542E+02 - 
.6926+01 - 
.2326+02 - 
.1256+03 - 
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MIN 
.1556+01 
>.2806 + 01 
.6066+01 
. 5726+01 
,1056+02 
.1166+02 
.7646+01 
« 9666+01 
.126E+01 
.1166+02 


MIN 
.351E+02 
.494E+02 
.4106+02 
. 1166+03 
.1676+03 
.1266 +0 3 
. 8916+02 
.1086+03 
.6926+02 
.1676+03 


J. Drag Coefficient = .00015 (SP9) 

Reducing the drag coefficient by one order of magnitude 
does not appear to greatly decrease the quality of the fore- 
cast in 24 hours. Figure IV-30 shows only slightly larger 
areas of departure from the base line than was seen in the 
last case, but these departures are of opposite signs. At 
500 mb, Figure IV- 31 shows only a small area of departure 
from the baseline. The statistical error summary for SF9, 
Table IV-IO, shows a generally slight increase in error 
at the surface and almost no difference at 500 mb in 
comparison with the baseline. 

More striking illustrations of the effect of tampering 
with the drag coefficient are seen in Figures IV-32 and IV-33 , 
which are the differences between the SF7 and SF9 runs at the 
surface and 500 mb, respectively. Here we can see differences 
of more than 8 mb at the surface and 60M at 500 mb. 
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FIGURE IV-31: DII'FERENCE BETt'TEEN BASELINE AND SF9 , SOOxnb 
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FIGURE IV-32: DIFFERENCE BETWEEN SP7 AND SF9, SURFACE 
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TABLE IV-10: 


ERROR STATISTICAL SUMMARY, SF9 
(DRAG COEFFICIENT « .00015). 


SURFACE 


?0M 

TO 

PHS 

SD 

MEAN 

MAX MIN 

0 

10 

.710E+00 

. 675E+00 

. 221E+00 

,2436+01 -.1566+01 

10 

20 

.118E+01 

. 118F+01 

-.4X9E-01 

.4016+01 -.286E+01 

20 

30 

.290E+01 

. 282E+01 

.679E+00 

.9316+01 -.5656+01 

30 

40 

, 357E+01 

.3196+01 

.1586+01 

.1126+02 -.64CE+01 

40 

50 

.414E+01 

. 315E+01 

. 269E+01 

.138E+02 -.7476+01 

50 

60 

•455E+01 

.342E+01 

. 3006+01 

,135E+02 -.8396+01 

60 

70 

.335E+01 

.274E+01 

. 1936+01 

,7956+01 -.6866+01 

70 

30 

.364E+01 

.3446+01 

. 121E+01 

.867E+01 -.833E+01 

80 

90 

» 236E+C1 

.141E+01 

. 190E+01 

,5596+01 -.3026+00 

0 

90 

.259E+01 

,2446+01 

. 871E+00 

.136E+02 -.8396+01 


lOH 

0 

30 

20 

30 

40 

50 

60 

70 

80 

0 


TO 

RMS 

10 

. 1656+02 

20 

.266E+02 

30 

.2276+02 

40 

.3806+02 

50 

.4806+02 

60 

.336E+02 

70 

.251E+02 

80 

, 3886 + 02 

90 

.. 4336+02 

90 

.2866+02 


500 MB 


SO 

.147E+02 
.248E+02 
. 191E+02 
« 346E+02 
.471E+02 
. 33BE+02 
.238E+02 
.219E+02 
. 114E + 02 
.276E+02 


MEAN 

MAX 

MIN 

i 

i 

.7496+01 

.5426+02 

-.3546+02 

3 

1 

1 

.9706+01 

,5226+02 

-.4986+02 

. j 

. 1246+02 

.6356+02 

-.4366+02 

1 

.1566+02 

.1316+03 

-.1126+03 

: ■ 1 

. 8986+01 

.1156+03 

-.1606+03 


-.1566+00 

.8416+02 

-.1166+03 

■■ ] 

-.7946+01 

,5396+02 

6396+02 


-.3216+02 

.1116+02 

-.9976+02 

] 

i 

-.4186+02 

-.2406+02 

-.6286+02 

1 

.7376+01 

.1316+03 

-.1606+03 
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K. Flat Earth (SF13) 

The elimination of terrain was detrimental principally 
in the lov/er levels. The two basic types of errors intro- 
duced at the surface, seen in Figure IV-34, are: 

The unimpeded movement of pressure systems demon- 
strated by the low over southern Afghanistan moving to the 
northwest and the central Asia high spreading in all direc- 
tions. Those show up in both the surface error chart. 
Figure IV-36, as a series of errors through central Asia 
and on the baseline departure chart, Figure IV-38, in the 
same region. 

The lack of eddy low formations, as seen on the 
surface chart, off the east coast of the U.S. and east of 
Greenland. 

The effects of the flat earth at 500 mb are generally 
a reflection of deviations at lower levels. SF13 at 500 mb. 
Figure IV-35, bears a close resemblance to the baseline. 

The error chart. Figure IV-37, is about the same as the 
baseline except for the China-USSR border area where the 
errors are smaller than those of the baseline. The amount 
of difference between the baseline and SF13, Figure IV-39, 
is about 60M in central Asia. 
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The improved forecast at 500 mb due to the flat earth 
model is undoubtedly a transitory feature. As the forecast 
advances beyond 24 hours, the errors at the surface will 
grow more rapidly and will, in turn, be reflected by in- 
creasing errors in the upper levels. 

The statistical error summary for SF13, Table IV-11, 
shows an overall increase in errors at the surface and some 
decrease of errors at 500 mb. 
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FIGURE IV- 34: FORECAST, SURFACE, SF13 
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FIGURE IV-35: FORECAST, SOOmb, SF13 
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FIGURE IV-36: ERROR PATTERN, SURFACE, SF13 
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FIGURE IV-38: DIFFERENC 
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FIGURE IV-39: DIFFERENCE BETWEEN BASELINE AND SF13, SOOmb 
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TABLE IV-11: ERROR STATISTICAL SUMMARY, SF13 

(NO TERRAIN) . 





SURFACE 



FROM 

TO 

RMS 

SO 

MEAN 

MAX 


0 

10 

.716E+00 

.6746+00 

. 243E+00 

.242E+01 

— 

10 

20 

,127r40l 

.127E+01 

.1086+00 

,4066+01 

- 

20 

30 

. 332E+01 

, 308E+01 

.122E+01 

.lOOE+02 


30 

40 

.391E+01 

• 365 E+01 

,1416+01 

.1356+02 


40 

50 

.459E401 

.42OE+01 

.184E+01 

.1456+02 

— 

50 

60 

.435E+01 

. 324E+01 

.290E+01 

,1186+02 


60 

70 

.317E+01 

.280E+01 

.1476+01 

,6876+01 

- 

70 

8C 

. 36SE+01 

.356E+01 

. eilE + 00 

.102E+02 

— 

80 

90 

,395E+01 

.235E+01 

.3176+01 

.1086+02 


0 

90 

.278E+01 

.264E+01 

.877E+00 

.1456+02 





500 

MB 



FROM 

T3 

RMS 

SO 

MEAN 

MAX 


0 

in 

. 165E+C2 

. 146E+02 

.772E+01 

.5436+02 

— 

10 

20 

f 265E + 02 

. 245E+02 

.lOOE+02 

,5896+02 


20 

30 

. 221E+02 

.190E+02 

. 112E+02 

.7926+02 


30 

40 

• 336E+02 

. 324E + 02 

. 898E+01 

,1056+03 


40 

50 

i445E+02 

.442E+02 

493E+01 

,6006+02 


5 0 

60 

.344E+02 

.327E+02 

-.1066+02 

,6026+02 


60 

70 

.318F+02 

.273E+02 

163E+02 

.4776+02 


70 

PO 

.423E+02 

.218E+02 

— . 3636+02 

-.3176+01 

— 

to 

90 

. 390E+02 

.151E+02 

-.3596+02 

-.4756+01 


f» 

Si 

90 

.278E+02 

.274E+02 

.4386+01 

.1056+03 

*“ ( 

original page t. 





POOR 

— 

quality 






MIN; 
.IS'tE + Ol 
.276E+01 
, 610E+01 
,959£+01 
.895E+31 
.689E+01 
.7B0E+01 
.790E+0X 
• 869E + 00 
,959E+01 


HIN 
.351E+02 
.496E+02 
. 405E + 02 
, 116E + 03 
. 165E+03 ; 
.127E+03 I 
► 901E + 02 : 
.108E + 03 . 
.618E+02 
.165E + 03 ' 
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L. No Precipitation (SP14) 

The effects of "turning off" precipitation on the model 
are almost all in the lowest levels during the first 24 
hours. This generally results in the reduction in the 
intensification of lows. In Figure IV-40, almost every low 
at the surface is less intense than in the corresponding 
baselxne situation {see SF17 [section IV-N] , following, for 
a comparison with the low relative humidity situation) . 
Although this is a problem in some areas, such as the U.S. 
East coast low development, it tended to reduce the surface 
errors (Figure IV-41) as compared to the baseline. This 
reduction occurred because much of the baseline error was 
due to overly rapid system development and movement. By 
elimination of the precipitation mechanism, the energy of 
these systems was reduced, and less development and movement 
took place. The difference between the baseline and SF14 at 
the surface is seen in Figure IV-42. 

Table IV-12 shows, inter alia, that for the entire 
northern hemisphere the RMS error is 2.42 mb for the surface 
SF14 vice 2.5 mb for the baseline, and that for SF14, the 
absolute maximum error shifted from positive to negative. 

All of the above does not prove the necessity for 
eliminating precipitation from the model. The apparent 
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salubrious effects discussed above are fleeting, and as the 
model progresses through time, the lack of the precipitation 
mechanisms would lead to gross underdevelopment of systems. 









FIGURE IV-41: ERROR PATTERN, SURFACE, SF14 
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FIGURE IV-42: DIFFERENCE BETWEEN BASELINE AND SF14 , SURFACE 
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TABLE IV-12: ERROR STATISTICAL SUMMARY, SF14 

{NO PRECIPITATION) . 


SURFACE 


FSOH 

TU 

PMS 

SD 

MEAN 

MAX 

MIN 

0 

10 

.718E+00 

.674E+00 

.248E+C0 

.247E+C1 

152E+01 

10 

20 

.118E+01 

.116E+01 

.220E-01 

.411E+C1 

282E+01 

20 

30 

.276G+01 

.272E+01 

.467E+00 

. 779E+01 

775E+01 

30 

^0 

.3HE+01 

.306E+01 

. 543E+U0 

. lOlE+02 

760E+01 

^0 

50 

.357E+01 

. 343E+01 

.999E+C0 

. 8E9E+01 

-.125E+02 

50 

60 

.^32E+01 

. 383E+01 

.203E+C1 

•102E+02 

131E + 02 

60 

70 

.351E+01 

.2eeE+oi 

.200E+01 

.601E+01 

717E+01 

70 

80 

.390E+01 

.370E+01 

.123E+01 

.933E+01 

-,937E + f>l 

80 

90 

.262E+01 

.159E+01 

.208E+01 

.609E+01 

780E+00 

0 

90 

. 2^2E+Q1 

. 236E+01 

. 545E+00 

,102E+02 

-.131E+02 





500 

MB 


FROM 

ro 

RMS 

SO 

MEAN 

MAX MIN 

0 

10 

.166E+02 

. 147E+02 

. 774E-I-01 

.551E+02 -,348E+02 

10 

20 

.270E+02 

«249£+02 

. 1C5E402 

.629E+02 -♦492E+02 

20 

30 

. 242E+02 

. 192E+02 

. 148E + 02 

.t44E+02 -,412E+n2 

30 

40 

.373E+02 

.329E+02 

. 176E + C2 

.129F+03 -.1O9E+03 

^0 

5 0 

.470F+02 

.461E+02 

.908E+01 

.llOE+03 -,155E+03 

50 

60 

.329E+02 

.327E+02 

. 362F + 01 

.736E+0H -.1C6E+03 

60 

70 

,244E+02 

. 240E+02 

-.403E+01 

.578E+02 -,8?2£+02 

70 

80 

.367E402 

.224E402 

290E + 02 

.125E+02 -.975E+02 

80 

90 

.406E+02 

.117E+02 

389E+02 

-.201E+02 -.621E+02 

0 

90 

.205E+O2 

. 271E+02 

. 873E+01 

.129E+03 -.155E+03 
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M. Balance Wind Initialization (SF16) 

The initialization of the winds by solution of the 
balance equation rather than by the analysis of wind fields 
produced a meteorologically good forecast that resembles 
the baseline in most respects. The one notable exception at 
the surface was in the development and movement of the East 
Asia low, as seen in Figure IV-43. SF16 avoided the false 
development of the trough west of Kamchatka and thereby 
eliminated the error in that area. It failed however, to 
move the main center of the low sufficiently eastward which 
lead to a new area of error as seen in Figure IV-44. Figure 
IV-45 shows that there are only slight differences between 
the baseline and SF16 elsewhere at the surface. 

At 500 mb. Figure IV-46, the departures from the base- 
line are even more subtle and tend to be associated with 
the East Asia low position and intensity as well as the 
slope of the trough south of Korea. The error pattern. 
Figure IV-47 , shows some variations from the baseline in 
that the errors are less severe over China and slightly 
greater south of Kamchatka and over Korea. This variation 
shows up on Figure IV-48. 

Although the summary of errors for SF16 (Table IV-13) 
shows a greater degree of error than the baseline, the 
difference is not significant. 


iy-71 



iy-72 


















CQ^GIPISIJ KSGE U 
OF POOR QUALTT^' 


3;y-74 








rV-76 












TABLE IV-13: ERROR STATISTICAL SUMMARY, SP16 

(BALANCED WINDS) . 


SURFACE 


FROM 

TO 

RMS 

SD 

MEAN 

MAX MIN 

0 

10 

.7D1E+00 

. 635E+00 

, 297E + 00 

,256E + 01 -,121E-f01 

10 

20 

.123E+01 

•iiie+oi 

, 520E + 00 

.3A7E+01 -.2A3E+CI 

zo 

30 

.259E+01 

.232E+01 

. 115E + Q1 

,013E+Ol -.A96E+01 

30 

AC 

.3A2B+01 

.305E+03 

.15AE+01 

.109E+02 -.616E+01 

AO 

50 

.A33E+01 

.399E+01 

. 168E + 01 

.127E+02 -.105E+02 

50 

60 

.A13E+01 

.310E+01 

.273E+01 

.112E+02 -.55AE+01 

60 

70 

.367E+01 

.29AE+01 

. 220F + 01 

»952t+01 -.70AE+01 

70 

80 

. A12E+01 

* A02E+01 

.901E+00 

.lllE+n2 -,863E+01 

80 

90 

.2A8C+01 

.21AE+01 

,126Ef01 

.551E+01 -.321E+01 

0 

90 

.25AE+01 

.235E+01 

.982E+00 

.127E+02 -.105E+02 





500 

MB 


FROM 

TO 

RMS 

SO 

MEAN 

MAX MIN 

u 

10 

.169E+02 

. 1A0E+O2 

. 815E+01 

.60AE+02 -.263E+02 

ID 

20 

•305E+02 

. 25AE+02 

. 169E+02 

•738E+02 -,380E+02 

2'J 

30 

.295E+D2 

. 213E+02 

.20AE + ')2 

.812E+02 -,AlIE+02 

30 

AO 

.339E+02 

. 311E+02 

. 13AE+02 

.955E+02 -.903E+02 

AO 

50 

.508E+P2 

, 507E + C2 

. 263E+U1 

.lGlE+03 ~.17Ab+03 

50 

60 

.362E+02 

. 362E+02 

. 23eE+Gl 

.691E+02 -.130E+03 

60 

70 

.255E+02 

.2A7E+02 “ 

• 606E+C1 

.A09E+Q2 -.768E+02 

70 

80 

.A38E+02 

.258E+02 - 

. 253E+02 

.1296+02 -.lOAE+03 

eo 

90 

.A82E+02 

.168E+02 - 

. A56E+02 

-.9A9E+01 -.776E+02 

0 

90 

.307E+02 

.291E+02 

.968E+01 

.lOlE+03 -.17AE+03 
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N. Relative Humidity at 25% (SP17) 

Rather than determine the initial relative humidity on 
the basis of system intensity. SF17 set the relative humidity 
at 25% over the entire map at the three lowest levels. 

The effect of this initialization was noticed primarily 
at the surface since the normal initial condition at the 
highest moisture level is of about the 25% magnitude. If 
the forecast extended much beyond 24 hours, however, changes 
in the upper levels would have become evident. 

The majority of the errors in the baseline were due to 
a combination of incorrect central pressures and system 
locations. SF17 did not greatly change system location but 
had an inhibiting effect upon developing systems, as can be 
seen in Figure IV-49 where the central pressures of the Gulf 
of Ala-ka low and the East Asia low are both about 4 mb higher 
than in the baseline. This can also be seen in Figure IV-50 
by the reduced error in these locations, and in Figure IV— 51 
where the difference between the baseline and SP17 is at least 
4 mb in several locations. These reductions in error are due 
principally to a loosening of gradient rather than better 
forecasts. 

The differences betv/een the baseline and SF17 at 500 mb 
are generally less than 30M. The statistical error summary for 
SF17, Table IV-14, shows a slight reduction in the overall RI4S 
and a shift of the absolute maximum to the negative. 
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TABLE IV-14: ERROR STATISTICAL SUMMARY, 

(RELATIVE HUMIDITY = 25 %) . 


SURFACE 


FROM 

TQ 

RKS 

SD 

0 

10 

.717E+00 

. 673E+00 

10 

20 

. li0E+Ol 

.118E+01 

^0 

30 

.278E+01 

, 274E+01 

30 

40 

.309E+01 

.303E+01 

40 

50 

,359E+01 

. 343E+D1 

50 

60 

. 430E+01 

. 378E+01 

60 

70 

. 350E+01 

.2F-7E+01 

70 

80 

,388E+01 

.369F+01 

fO 

90 

.262E+01 

.159E+01 

0 

90 

.242E+D1 

. 236E+01 


500 


FROM 

TO 

RMS 

SD 

0 

10 

.166E+02 

. 147E+02 

10 

20 

.270E+02 

.2498+02 

20 

30 

.243E+02 

. 193E+02 

30 

40 

.373E+02 

.329E+02 

40 

50 

.471E+02 

.463E+02 

50 

60 

.328E+02 

• 3 26 E +0 2 

60 

70 

.245E+02 

. 241F. + 02 

70 

80 

.367E+02 

. 224E+02 

60 

90 

.408E+02 

.118E+32 

0 

90 

• 285E+02 

. 272E+02 


MEAN 

MAX 

MIN 

.247E+00 

.246E+01 

-.153E+01 

.155E-01 

.409E+01 

-.285E+01 

. 486E+00 

.793E+D1 

812E + 01 

. 585E+00 

.102E+02 

721E+01 

.104E+01 

.874E+01 

123E+02 

.205E+01 

.105E+02 

-.127E+02 

.201E+01 

.799E+01 

712E+01 

. 121l; + 01 

.934E+01 

-.940E+01 

. 208E+01 

. 606 E + 01 

855E+00 

.555E+00 

.105E+02 

127E+02 

MB 

MEAN 

MAX 

HIN 

.773E+91 

.551E+02 

-.346E+92 

. 105Ef02 

.630E+02 

-.492E+C2 

. 148E+02 

.649E+C2 

414E+C2 

• 176E + 02 

.128E+03 

-.109E+03 

.697E+Q1 

.llOE+03 

-.1566+03 

. 346E + 01 

.741F+02 

108E+03 

. 419E+01 

,583E+02 

-.823E+02 

.291E+02 

.125F+02 

-.973E+02 

. 390E + 02 

-,20QE+02 

-.622E+02 1 

. 868E + 01 

.128E+03 

156E+03 


EPOR 
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O. Relative Hiimidity at 75% (SF18) 

The over-specification of initial relative humidity in 
SF18 introduced large errors in two surface locations (see 
Figure IV-52) . The East Asian low, which is undergoing rapid 
development and movement, has these processes accelerated. 

Where there is some overdevelopment of troughs to the east 
and north of this system in the baseline, SF18 has a 975 mb 
low form in the east and a 1000 mb low form to the north. 

The other location is over Worth Africa. What was merely a 
trough over Tunisia on the baseline became a 1000-mb low 
when excessive moisture was introduced. 

Figure IV-53 shows surface error patterns for SF18. 

Most of the errors are greater than for the baseline, but a 
few selected systems, such as the Gulf of Alaska low, show a 
decrease in error. The difference between the baseline and 
SP18 at the surface. Figure IV-54, clearly shows the selectivity 
of the departure from baseline. The differences between the 
baseline and SF18 at 500 mb were generally less than 30M. 

As would be expected, the error suramary, Table IV-15, 
shows fairly sizeable increases in error at the surface. 

One should note that overspecifying the amount of 
moisture, especially in rapidly developing systems, can 
greatly affect the model forecast. 
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TABLE IV- 15 


: ERROR STATISTICAL SUMMARY, SF18 

(RELATIVE HUMIDITY = 75%) . 

SURFACE 


FRDH TO 
0 10 
10 20 
20 30 

30 ^0 

tfO 50 
5C 60 
60 70 
70 BO 
8C 90 
0 90 


FROM TO 
0 10 
10 20 
20 30 

30 1*0 

40 5 0 

50 60 

60 70 

70 00 

BO 90 
0 90 


RMS 
.706E+00 
. 116F+01 
.294E+01 
• 406E+01 
.410E+01 
.440E+01 
,322E+01 
, 369E+01 
. 227E+01 
.265E+01 


RMS 
.164E+02 
.255E+02 
.22‘6E + 02 
♦364E+02 
.477E+02 
. 347E+02 
.257E+02 
.402E+02 
.441E+72 
.284E+Q2 


SD 

,673E + nr) 
,il6E+01 
. 269E + 01 
.333E+01 
. 324E+01 
. 365E+01 
t 27 6E+P1 
.360E+01 
.155E+01 
.247E+01 


500 

SD 

.147E+02 
,249E+02 
. 196E+02 
. 345E+02 
.476E+02 
,346E + ''2 
.238E+02 
.222E+C2 
. 115E+02 
.278E+02 


MEAN 
.214E+00 
,216F-01 
. 117E+01 
.232E+01 
.25 IE + 01 
, 246E+01 
.167E+01 
.783E+00 
. 166E+01 
.972E+00 


MB 

MEAN 
.741E+01 
.915E+01 
. 1136+02 
. 114E+02 
. 325E+01 
334E+01 
963E+C1 
335E+02 
426C+02 
. 580E+01 


MAX 
•241E+01 
.423E+01 
.983E+01 
.146E+02 
.972E+01 
.170E+02 
.73GE+C1 
, 8a4E+01 
.557E+01 
.170E+C2 


MAX 
.538E+02 
,625E+02 
.636E+02 
, 117E + 03 
. 109E + 03 
. 109E+03 
.534E+02 



MIN 

-,162E 

+ 01 

“,292E 

+01 

-.483E 

+ 01 

-.578E 

+ 01 

-.935E 

+ 01 

-.941E 

+ 01 

-.721E 

+ 01 

-,924E 

+01 

-. 119E 

+ ul 

-.941E 

+ 01 


MIN 
356E+02 
501E+02 
487E+02 
112E+03 
164E+03 
119E+03 
6O9E+02 
E + 03 
E + 02 
E+03 


♦911E+01 -.102 
-.242E+02 -.650 
.U7E + 03 -.164 
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P. Modified Truncation/No Terrain Tapering (SP19) 

SF19 was run with two computational devices modified. 

One modification was to reduce the truncation tendency linearly 
from 1 to . 5 between 20°N and 5°N and then have a constant .5 
south of 5°N. The other modification is to retain the ful?. 
(smoothed) terrain heights up to the map borders. 

Not surprisingly, SP19 produced large areas of undesir- 
able patterns. Figure IV-55 shows that at the surface these 
are evident south of 20®N. Figure IV-56 shows that these 
patterns introduced errors as large as 12 mb. 

In spite of the fact that the devices were only modified 
south of 20°N, the departure from the baseline extended north 
of 60 °N in 24 hours as is evident in Figure IV-57. In 72 
hours, the entire domain would be strongly disturbed by these 
modifications . 

At 500 mb. Figure IV-58, some irregularities are observed 
below 20°N that do not appear to be as intense as those seen 
at the surface. Figure IV-59 illustrates that these errors 
are not much more than 60M greater than the errors of the 
baseline. This slight departure is even more readily apparent 
in Figure IV-60. 

Although both truncation and terrain modifications were 
introduced to SF19, it appears that the truncation factor was 
far more critical. This is assumed by the facts that: terrain 
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does not influence greatly the upper levels but in SF19 the 
500 mb is affected; the area below 20°N is predominantly 
sea level; and areas over land do not appear to be more 
adversely affected than do the ocean areas. 

In Table IV-16 , the Rt4S errors are greater than the 
baseline at the surface at every latitude band. At 500 mb, 
it exceeds the baseline up to 60°N. 
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FIGURE IV-58: FORECAST, 500mb 


IV-J4 












: I 


. \ 


- page is 

DE POOR* Q UACma 


IV- 9 6 





TABLE IV-16; ERROR STATISTICAL SUMIIARY, SF19 
{MODIFIED TRUNCATION/TERRAIN 
TAPERING) . 


SURFACE 


FROM 

TO 

RMS 

SO 

MEAN 

MAX 

HIM 

0 

10 

.430E+O1 

.43AE+01 

-.524E+CO 

.954E+01 

838E+01 

10 

20 

.332E+01 

. 328E+01 

546E + 00 

.658E+01 

-.653E+01 

20 

30 

.366E+01 

.363E+01 

. 484E+00 

, 799E+01 

-.755E+01 

30 

40 

.373E+01 

. 36PE+01 

, 980E+00 

.106E+02 

-.615E+01 

40 

50 

.335E+01 

.346E+D1 

, 169E+01 

.llXE+02 

-.914E+01 

50 

60 

.461E+01 

.374E+01 

.271E+01 

.il5E+02 

112E+0Z 

60 

70 

.336E+01 

. 254E+01 

.220E+01 

.S59E+01 

-.641E+01 

70 

30 

, 376E+01 

.350E+01 

. 137E+01 

.982E+01 

823E+01 

eo 

90 

.241E+01 

.149E+01 

. 189E + 91 

,589E+01 

-*.661E+00 

0 

90 

.389E+01 

.387E+01 

. 355E+00 

.115E+02 

112E+02 





500 

MB 



FROM 

TO 

RMS 

SD 

MEAN 

MAX 

MIN 

3 

10 

*4628+02 

.460E+02 

.4346+01 

.953E+02 

-.7636+02 

10 

20 

.421E+32 

. 415E+02 

. 751E+01 

.9146+02 

-.6336+02 

20 

30 

.307E+02 

,2956+02 

. 829E+01 

.7196+02 

6636+02 

30 

40 

.3978+02 

. 389E+02 

. 791E+0I 

.1206+03 

-.1156+03 

40 

50 

.4956+02 

.4956+02 

, 120E+C1 

.1106+03 

-.1636+03 

50 

60 

.378E+02 

. 378E+C2 

-.187E+ni 

.691E+02 

-.1196+03 

60 

70 

.241E+02 

.234E+02 

571E+01 

.603E+02 

-.0126+02 

70 

80 

,3956+02 

.2266+02 

-.3106+02 

,642E+01 

-.9766+02 

80 

90 

.437E+02 

.108E+02 

-.424E+02 

-.2656+02 

6106+02 

0 

90 

. 414E+02 

.412E+02 

.3826+01 

.1206+03 

-.1636+03 
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Q. Multiple Modification (SFll) 


The previous studies generally resulted in subtle 
changes in the forecast. In addition to those, it was 
instructive to demonstrate the effect of the modification of 
several devices simultaneously. In SFll, the coefficients 
for temperature and momentum diffusion, as well as for surface 
drag and pressure smoothing, were set to zero. In addition, 
the precipitation and heating processes were eliminated. This 
would be an adiabatic, frictionless version of the standard 
model. 

In 24 hours, the damage to the forecast was drastic. At 
the surface, Figure IV-61, most of the systems are recogniz- 
able but all have undergone excessive development. The 
pattern of Figure IV-62 has errors in excess of 20 mb and 
almost without exception there are large departures from the 
baseline in Figure IV-63. 

At 500 mb. Figure IV-64, systems have assumed non- 
meteorological gradients and shapes. While some of the errors 
in Figure IV-65 resemble those of the baseline, there are new 
regions with large errors. Figure IV-66 shows the departure 
from the baseline to be largely confined to Asia although 
other departures exist. 

Table IV-17 indicates large errors associated with SFll. 
The maximum errors were 40 mb at the surface and 414M at 

500 mb, while the baseline was 12.7 mb and 163M at those 
levels. 
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FIGURE IV-6I: FORECAST, SURFACE, SPll 


Qt POOR 


IV- 9 9 
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TABLE IV-17; ERROR STATISTICAL SUMMARY, SFll 
(MULTIPLE MODIFICATIONS) . 


SURFACE 


FRQH 

TO 

RMS 

SD 

MEAN 

MAX 

MIN 


0 

10 

.727E+00 

, 700E+0Q 

. 196E + 00 

.256E+01 

“.152E+01 


10 

20 

. lAAF+01 

. lAAE+01 

-» 106E + 00 

.A69E+0X 

AOAE + 01 


20 

30 

,387E+01 

-.381E + 0 J 

. 657E+00 

.1296+02 

1A1E + G2 


30 

AO 

.63AE+01 

.621E+01 

.125E+01 

.AOOE+02 

-» 176E+02 


AO 

50 

466 AE+OI 

*595E+01 

. 295E + C 1 

.225E+02 

151E + 02 


50 

60 

. 625E+01 

. 552E+0L 

.292E+01 

.255E+02 

105E+02 


60 

70 

.391E+01 

, 3A3E+Q1 

.183E+01 

.926E+01 

-.792F+01 


70 

60 

. AA9F+01 

, A36E+01 

« 105E+01 

.977E+01 

-.955E+D1 


• 80 

90 

.215E+01 

. 2 C 6 E +01 

* 6 A 2 E + 00 

.622E+01 

“» 386E + 01 

/ 

0 

90 

. 38AE+01 

.375E+01 

.611E+00 

. A00E+C2 

-.176E+02 

.] 




500 

IVIB 



i 

FRQH 

TD 

RMS 

SD 

MEAM 

MAX 

MIN 


0 

10 

. 152F+02 

.1A6E+02 

.599E+01 

.557E+02 

3A5E+02 


10 

20 

.279E+02 

. 266E+02 

.833E+01 

♦701E+02 

~»A06E+OH 


20 

30 

.355E+02 

.299E+02 

. 191E + 02 

.162E+03 

-.620E+02 


30 

AO 

.7DAE+02 

. 632E+02 

.308E+02 

. A lAE + 03 

115E+03 


AO 

50 

♦ 637E+P2 

.56AE+P2 

. 255E+02 

, 262E + 03 

1A6E+03 


50 

60 

.A27E+02 

^AilE +02 

. llAE->02 

.1706+03 

955E + 02 

1 

fcO 

70 

.288E+02 

.286E+02 

-.3A9E+01 

. 6 8 AE + 02 

875E + 02 

■ i 

70 

80 

.381E+02 

.255E+02 

263E+02 

. 361E + 02 

108E+03 

. f 

80 

90 

.337E+02 

, 186E+02 

“,2f lE+02 

-.53PE+00 

-.6A7E+02 


0 

90 

.392E+02 

. 372E+02 

. 122E+02 

. AlAE + 03 

1A6E + 03 
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R. No Heating (SF20) 


In the final run of the model, the processes that 
generate heating contributions every hour were eliminated. 
The departure of this run from the baseline was negligible. 
The greatest departure from the baseline was less than .5 
mb at the surface and 1.5M at 500 mbs. 


V. COMPARATIVE PERFORMANCE; ENERGETICS, PRECIPITATION 

AND STORM INTENSITY 

Characteristics, effects, devices, and/or coefficients 
peculiar to the standard forecast model (SP-3) v/ere altered, 
singly or in combination, in order to determine the relative 
impacts on a one-day forecast. Once presented, these impacts 
can be ordered and discussed. 

A. Energetics 

The effects of each model alteration on 24-hour forecasts 
of layer-mean kinetic energy and square vorticity are tabulated 
and described. 


1. Temporal Filtering 

Tables V— 1 and V— 2 show the effect of temporal 
filtering on forecast changes of mean kinetic energy and mean 
square vorticity, respectively. In both, the impact is shown 
to be negligible. 

2. Radiation 

Tables V— 3 and V— 4 show the effect of radiation cal- 
culations on forecast changes of mean kinetic energy and mean 
square vorticity, respectively. In both, the impact is negli- 
gible. 


3. 


Vertical Resoln^-i nn 


Tables V-5 and V-6 show the effect of doubling the 
vertical resolution (from five to ten layers) on forecast 
changes of mean kinetic energy and mean square vorticity, 
respectively, m both cases, doubling the resolution leads to 

greater losses in low levels and to greater gains at high 
levels. 


Temperature Diffusion 

Tables V-7 and V-8 show the effect of temperature 
on forecast changes of mean kinetic energy and mean 
square vorticity, respectively. The coefficient (TDK) was 
reduced by an order of magnitude (from the standard model) . 
This led to increases in forecast KE of about 5%, generally. 
The impact on square vorticity varied with level, with maximum 
impact near the Sigma=0.3 level (strongest winds). 

^ • Init:xa.l Winds 

Tables v-9 and V-10 show the effect of a variation 
in the initial wind specification on forecast changes of mean 
kinetic energy and mean square vorticity, respectively. The 
atandard model used analyzed winds. Run SF-16 was initialized 


with non divergent winds obtained through solution of the 
non-linear balance equation. In general, the balanced winds 

produced smaller increases (or larger decreases) in predicted 
KE and SV. 

6 . Precipitation 

Tables V-11 and V-12 show the effect of precipitation 
on forecast changes on mean kinetic energy and mean square vor- 
ticity, respectively. As expected, smaller increases (or 
larger decreases) in predicted KE and SV are indicated in the 
moisture-bearing levels • 

Initial Moisture Specification 

Tables V-13 and V-14 show the effect of variations 
in the initial moisture specification on forecast changes of 
mean kinetic energy and mean square vorticity, respectively. 
Examine Runs SP-17 and SP-18, first. There is a difference of 
50-O in the initial relative humidities. The greater RH pro- 
duced larger increases (smaller decreases) in KE and SV, 
generally. The effect seems to attenuate with altitude (level) , 
but there are exceptions. The 28% difference in SV at Sigma=0.3 
level is anomalous. Next, note that the actual run (SP-3) seems 
to indicate that its (initial) humidity condition v/as bounded by 
these extreme conditions, mainly. 


V-3 


8 . 


Mountains 


Tables V-15 and V-16 show the effect of mountains 
on forecast changes iii mean kinetic energy and mean square 
vorticity, respectively. Clearly, mountains lead to larger 
increases in both KE and SV, mainly as a consequence of in- 
creased latent heating. The effect seems to attenuate with 
height through the moisture-bearing layers. Minor alterations 
to the mountains (Run SF-19) in the low latitudes had minor 
effect (there is precipitation south of about 25°N in these 
models) . 


9 . Pressure Smoothing 

Tables V-17 and V-18 show the effect of variation in 
the non-linear pres sura- smoothing coefficient (SM) on forecast 
changes of mean kinetic energy and mean square vorticity, res- 
pectively. This effect is rather large, and it decreases with 
height. The KE and SV forecast changes are larger positive 
(smaller negative) as SM is reduced in magnitude. 

10. Momentum Diffusion 

Tables V-19 and V-20 show the effect of momentum 
diffusion on forecast changes in mean kinetic energy and mean 
square vorticity, respectively. In Run SF-2, the coefficient 


V-4 


(DK) was reduced by an order of magnitude. This led to large 
increases in the forecast changes in KE and SV. Indeed, when 
the term is caused to vanish (Run SF-10) , extremely large 
increases result. The term is obviously needed to maintain 
reasonable levels of KE and SV, but not so strong that net 
decreases result (Run SP-1) . 


11. Surface Friction 


Tables V-21 and V~22 show the effect of friction on 
forecast changes in mean kinetic energy and mean square vor- 
ticity, respectively. Note that the tests were conducted for 
two magnitudes of pressure smoothing. 

In general, the two-fold increase in friction led to 
smaller increa.^es (larger decreases) in forecast KE and sv. 

At sigma=0.9, the KE difference is 14.6 when SM=0 . 01 and 11.1 
when SM-0.02. At Sigma=0.9, the SV difference is 19.3 when 
SM=0.01 and 14.4 when SM=0.02. Thus, the 100% variation in 
riction produces different results when an important dis- 
cretionary procedure such as pressure smoothing is operative. 


B. 


Storm Intensity 


Eight extratropical low-pressure centers V'/ere present 
at 1200Z, 22 April 1976. The impact of each discretionary 
procedure/effect tested on the predicted intensity of these 
pressure systems is shown in Table V-23. 

By far, the 3X variation in horizontal grid resolution 
had the greatest impact (probably because the precipitation 
doubled as the resolution tripled) . Precipitation (3MB) is 
also important, but momentum diffusion can be equally important 
to the outcome. Also, pressure smoothing (2.5MB) has equal 
impact to a 50% variation in the initial moisture specifi- 
cation (2.6MB). Both radiation and temporal filtering had 
no effect on storm intensity. 

C. Precipitation 

The predicted amounts of precipitation (both large-scale 
and convective-scale) north of 25°N were tabulated for each 
model run. 

Large variations ('^100%) occurred as a consequence of 
tripling the horizontal resolution or of a 50% variation in 
the initial moisture specification. Doubling the vertical reso- 
lution from five to ten layers resulted in a 10-15 percent varia- 
tion in precipitation only. All other parameters/effects tested 
led to only minor variations in precipitation (less than five 
percent, generally) . 


V-6 


TABLE V-1: 


EFFECT OF TEMPORAL FILTERING ON FORECAST 
CHANGES OF MEAN KINETIC ENERGY (IN 
PERCENT) . 


RUN 



SIGMA 

LEVEL 


0.9 

0.7 

0.5 

0.3 

0.1 

FILTER 
ACTIVE 
(RUN SF-3) 

18.1 

1.0 

1.9 


-2.6 

\ 

-23.3 

FILTER 
INACTIVE 
(RUN SF-15) 

18.2 

1.0 

1.9 


-2.6 

j 

i 

i 

i 1 

i -23.3 ^ 

f j 


NOTE: 


THESE ARE 2 4 -HOUR FORECAST CHANGES 







TABLE V-2; 


EFFECT OF TEMPORAL FILTERING ON FORECAST 
CHANGES OF MEAN SQUARE VORTICITY (IN 
PERCENT) . 


RUN 

'r‘ "■ — —i 


SIGMA 

LEVEL 



1 

0.7 

0.5 

i 0.3 

0.1 i 

1 

FILTER 
1 ACTIVE 
i (RUN SF-3) 

1 ) 
i 

24.1 

-13.5 

8.0 

23.2 

j 

r 

-43.1 ; 

1 

FILTER i 

INACTIVE 
(RUN SF-15) 

1 

i 

24.3 i 

-13.5 

8.0 

i i 

i * 

i 23.4 , 

i 1 

-43.1 


NOTE: THESE ARE 24 -HOUR FORECAST CHANGES. 
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TABLE V-3: 


EFFECT OF RADIATION AND SENSIBLE HEATING 
ON FORECAST CHANGES OF MEAN KINETIC 
ENERGY (IN PERCENT) . 


CONDITION/RUN 

SIGMA LEVEL 

0.9 

0.7 

0.5 

0.3 

0.1 

NO RADIATION OR 
SENSIBLE HEATING 
(RUN SF-20) 

18.1 

0.9 

1.9 

-2.8 

-23.3 

i 

i 

STANDARD 
(RUN SF-3) 

i 

1 18.1 
j 

1.0 ' 

! 1.9 

i 

-2.6 

-23.3 

t 

1 








TABLE V^4: 


radiation and sensible HEATINr 

ON FORECAST CHANGES OF MEAN SOUARF 
VORTICITY (IN PERCENT) . ^ 


CONDITION/RUN 


SIGf4A LEVEL 


NO RA 'STATION OR 
SENSIBLE HEATING 
(RUN SF-20) 


STANDARD 
(RUN SF-3) 


24.1 


24.1 


-13.5 


-13.5 


0.3 ! 


23.5 


23.2 


V-10 


TABLE V-5: EFFECT OF VERTICAL RESOLUTION ON FORECAST 

CHANGES OF MEAN KINETIC ENERGY (IN 
PERCENT) . 


CONDITION/RUN 


SlGfJlA LEVEL 


5-LEVELS 
(RUN F18) 


“4.9 


-13.1 j -10.9 -19.5 


10-LEVELS 
(RUN T4) 


"12.5 -13.8 


-18.4 


DIFFERENCE 


-3.5 


-2.9 


WITH A DIFFERENT 


those used for the other sensitivity*^ 


TESTS 


TABLE V-6: EFFECT OP VERTICAL RESOLUTION ON FORECAST 

CHANGES OP MEAT^ SQUARE VORTICITY (IN 
PERCENT) . 


CONDITION/RUN 

SIGMA LEVEL 

0.9 

0.7 

0.5 

0.3 

0.1 

5 -LEVELS 

-12.2 

! 

-43.2 

-37.1 

-37.4 

' --43.4 

(RUN P18) 





: 

1 

10 -LEVELS ! 

-25.3 

-44.2 

-44.3 

-33.5 

' 

-21.0 

(RUN T4) I 




1 

1 I 

! I 

1 

DIFFERENCE 

-13.1 

i 

-1.0 1 

-7.2 

i 

3.9 

" ■ ' ■' r 

22.4 


NOTE: THESE RUNS WERE INITIALIZED WITH A DIFFERENT SET OF 

ANALYSES THAN THOSE USED FOR THE OTHER SENSITIVITY 
TESTS. 
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TABLE V-8; EFFECT OF TEMPERATURE DIFFUSION ON 
FORECAST CHANGES OF MEAN SQUARE 
VORTICITY (IN PERCENT) . 


COEFFICIENT/RUN 



SIGMA 

LEVEL 

1 

0.9 

0.7 

0.5 

0.3 

0.1 i 

TDK = 10^ 

23.3 

-9.6 

24.0 

44.3 

-41.0 : 

(RUN SF-4) 





i 

TDK = 10^ 

i 

1 

24.1 

-13.5 

8.0 

23.2 

r 

-43.1 : 

] (RUN SF-3) 

1 ■ 



V 


j 

DIFFERENCE 

1.2 

3.9 

16.0 

21.1 

2,1 


V-14 





TABLE V-9: 


EFFECT OF INITIAL WIND SPECIFICATION ON 
FORECAST CHANGES OF I-lEAN KINETIC ENERGY 
(IN P.ERCENT) . 


CONDITION/RUN 

1 


SIGMA 

LEVEL 


0.9 

0.7 

0.5 

0.3 

0.1 

ANALYZED WINDS 
(RUN SF-3) 

: 18.1 

1.0 

1.9 

-2.6 

-23.3 

BALANCED* WINDS 
(RUN SF-16) 

1 

' 5.4 

I 

-0.4 

-3.1 

! 

-7.5 

-19.4 I 
i 

DIFFERENCE 

j-12.7 

-1.4 

-5.0 

-4.9 

3.9 


* FROM NON-LINEAR BALANCE EQUATION; NO SPECIFICATION OF 
DIVERGENT COMPONENT. 


V-15 



TABLE V-10: EFFECT OF INITIAL WIND SPECIFICATION 

ON FORECAST CHANGES OF MEAN SQUARE 
VORTICITY (IN PERCENT) . 


CONDITION/RUN 

SIGMA LEVEL 

i 

0,9 

0.7 

0.5 

0,3 

■ 0.1 i 

ANALYZED WINDS 
(RUN SP-3) 

24.1 

-13.5 

i 

8.0 

23.2 

i -43.1 

i 

j 

BALANCED* WINDS 
(RUN SF-16) 

— i 

12.6 

-16.8 

-16.2 1 

[ 

' -10.5 

! -36.7 ! 

1 
1 

difference j 

-11.5 

-3.3 

-24.2 

-33.7 

6.4 


:k 


PROM NON-LINEAR BALANCE EQUATION- 
DIVERGENT COMPONENT. 


NO SPECIFICATION OP 


V-16 



TABLE V-11: EFFECT OF PRECIPITATION ON FORECAST 

CHANGES OF MEAN KINETIC ENERGY (IN 
PERCENT) . 



BOTH CYCLONE-SCALE AND CONVECTIVE-SCALE . 


-i: 





TABLE V-12; EFFECT OF PRECIPITATION ON FORECAST 
CHANGES OF MEAN SQUARE VORTICITY 
(IN .PERCENT) . 



V-18 










TABLE V-14: EFFECT OF INITIAL MOISTURE (RHo) SPECIFI- 

CATION ON FORECAST CHANGES OF MEAN SQUARE 
VORTICITY (IN PERCENT) . 


CONDITION/RUN 



SIGMA 

LEVEL 


0.9 

0.7 

0.5 

0.3 

0.1 

RH„ = 25% 
(RUN SF-17) 

-3.6 

-23.1 

5.1 

18.1 

-42.7 

RHo = f(VORT)* 
(RUN SF-3) 

24.1 

-13,5 

i 

8.0 

23.2 

-43.1 

RHo = 75% 
(RUN SF-18) 

27.7 

1 

-13.5 

17.7 

46.1 

1 

1 

i 

-41.8 1 

» 

I 

DIFFERENCE 
(75% - 25%) 

1 

31.3 

9.6 

12.6 

28.0 

0.9 


* INITIAL RH PARAMETERIZED IN TERMS OF RELATIVE VORTICITY AT 
EACH PRESSURE LEVEL. 






TABLE V^15* EFFECT OF MOUNTAINS ON FORECAST CH^GES 

table V 15. kinetic energy (IN PERCENT) . 


CONDITION/RUN 


SIGMA LEVEL 


MODIFIED* 
NEAR BORDERS 
{RUN SF-3) 


18.1 


- 2.6 


UNMODIFIED 
(RUN SF-19) 


14.1 -0.9 


FLAT EARTH 
(RUN SF-13) 


i -3.0 -5.4 


-4.4 


DIFFERENCE 21.1 

(SF-3 VS SF-13) 


* FLAT EARTH SOUTH OF 5°N; UNMODIFIED NORTH OF 20 “Nj 
TAPERED BETWEEN 5°N AND 20°N. 


TABLE V-16: 


EFFECT OF MOUNTAINS ON FORECAST CHANGES 
OF MEAN SQUARE VORTICITY (IN PERCENT) . 


CONDITION/RUN 

SIGMA LEVEL 


0.9 

I 0.7 

1 0.5 

0.3 

; 0.1 1 

MODIFIED* 
NEAR BORDERS 
(RUN SF-3) 

\ 

t 

24.1 

1 

-13.5 

! 8.0 

! 

1 

23.2 

1 -43.1 1 

! i 

I 

j ! 

UNMODIFIED 
(RUI'J SF-19) 

24.1 

-14.4 


23.6 

1 i 

[ 1 

-45.2 

1 

1 

1 

FLAT EARTH ! 

(RUN SF-13) j 

! 

! 

-16.5 

i 

1 

1 

-31.2 

0.0 

11.1 

-45,6 

1 

1 

DIFFERENCE 
(SF-3 vs SF-13) 

40.6 

17.7 

8.0 

12.1 

i 

2.5 


•k 


FLAT EARTH SOUTH OF 5°N; UNMODIFIED 
TAPERED BETWEEN 5°N AND 20°N. 


NORTH OF 20“N; 


V-22 


TABLE V-17: 


EFFECT OF PRESSURE-SMOOTHING ON FORECAST 
CHANGES OF MEAN KINETIC ENERGY {IN 
PERCENT) . 


COEFFICIENT/RUN 

SIGMA LEVEL 

0.9 

0.7 

0.5 

0.3 

0.1 

SM = 0.0 
(RUN SF-5) 

31.6 

6.2 

3.0 

-1.7 

-23.3 

SM = 0.01 
(RUN SP-3) 

18.1 

1.0 

1.9 

i 

-2.6 

-23.3 

SM ^ 0.02 
(RUN SF-6) 

1 

3.5 

-3.8 

0.8 i 

1 

1 j 

1 -3.5 

1 1 

-23.9 

DIFFERENCE 
(SF-3 VS SF-6) 

14.6 

1 

4.8 

1 

1.1 1 

1 

: 0.9 

I 0.6 


NOTE: THESE ARE 24-HOUR FORECAST CHANGES. 


V-23 



TABLE V-18: 


EFFECT OF PRESSURE SMOOTHING ON FORECAST 
CHANGES OF MEAN SQUARE VORTICITY (IN 
PERCENT) . 


i 



SIGMA 

LEVEL 

! 

COEFFICIENT/RUN ; 

0.9 i 

t 

0.7 

0.5 

0.3 i 

0.1 

i SM = 0 . 0 
(RUN SF-5) 

i 

57.8 

' 

-2.9 

10.9 

27.1 

1 

-42.7 

1 

1 

1 

i SM = 0.01 

1 (RUN SF-3) 1 

i ' 

i I 

24.1 

-13.5 

i 

1 

8.0 

23.2 

j 

-43.1 

1 1 

i SM = 0.02 i 

(RUN SF-6) 1 

1 i 

0.0 

i 

i 

! -21.2 

5.1 

t 

20.2 

1 

-42.7 i 

i ! 

DIFFERENCE 
(SF-3 vs SF-6) 

24.1 

7.7 

2.9 

; 3.0 

-0.4 


NOTE: THESE ARE 24-HOUR FORECAST CHANGES. 


V-24 


TABLE V-19: EFFECT OF MOMENTUM DIFFUSION ON 

FORECAST CHANGES OF MEAN KINETIC 
ENERGY (IN PERCENT) . 


COEFFICIENT/RUN 1 



SIGMA 

LEVEL 


1 

0.9 

0.7 

0.5 

j 0.3 

0.1 


DK = 

0.0 

112.0 

83.1 

91.5 

i 74.1 ’ 

t 

0.6 

f 

(RUN 

SF-10) 




1 1 

1 1 

! j 



DK = 

10^ 

o 

* 

in 

CO 

54.0 

54.4 

i 1 

42.6 ' 

- 9.3 


1 (RUN 

[ 

i 

SF-2) 


1 

r 

r 


1 

i 

j DK = 

10 ^ 

18.1 

i 1.0 

! 1.9 

i 

j ^2.6 

-23.3 

f 

1 

1 

1 

! 

(RUN 

1 

SF-3) 


1 


; 


1 

i DK = 

f (K) * 

o' • ‘J 

I 

1 

; - 15.9 

1 - 14.8 

^ -18 • 8 

■ -26.4 


j (RUN 

SF-1) 


1 

1 

1 

i 

1 

I 



* DK 


1.0 

E06, 

K=5 

2.2 

E06, 

K=4 

2.2 

E06, 

K=3 

1.9 

E06, 

K=2 

1.6 

E06, 

K-1 




V-25 


TABLE V-20: 


effect of momentum diffusion on 
changes on MEAN SQUARE 
VORTICITY (IN PERCENT) . 


COEFFICIENT/RUN 


SIGMA LEVEL 


DK = 0.0 
(RUN SF-10) 


172.1 139.3 


251.9 321 


DK = 10^ 
(RUN SF-2) 


131.4 


96.3 151.7 200. 


DK = 10° 
(RUN SF-3) 


24.1 


-13.5 


3K = f (K) * 
(RUN SF-1) 


-9.8 J -41.6 


-35.9 


DK = 


1.0 E06, K=5 
2.2 E06, K=4 
2.2 E06, K=3 
1.9 E06, K=2 
1.6 E06, K-1 




TABLE V-21: 


EFFECT OF SURFACE FRICTION ON FORECAST 
CHANGES OF MEAN KINETIC ENERGY (IN 
PERCENT) , UNDER TWO CONDITIONS OF 
PRESSURE SMOOTHING. 


COEFFICIENT/RUN 



SIGMA 

LEVEL 

! 

0.9 

0.7 

1 0.5 

0.3 

0.1 

1 

1 PRESSURE-SMOOTHING 

1 

COEFFICIENT = 0.01 

I 

CD = 0.00015 
(RUN SF-9) 

38.0 

i 

] 

8.5 

5.8 

i 

I 

-0.7 

i 

-22.6 

j 

CD = 0.0015 
(RUN SF-3) 

i 

18.1 

1 

1.0 

1.9 

i 

-2.6 

I 

-23.3 

1 

CD = 0.003 
(RUN SF-7) 

3.5 

-5.7 

-1.7 

-4.7 

1 

j 

-24.5 

PRESSURE-SMOOTHING 

COEFFICIENT =0.02 

1 

CD = 0.0015 
(RUN SF-6) 

3.5 

-3.8 

0.8 

-3.5 

-23.9 

1 CD = 0.003 
(RUN SF-8) 

-7,6 

-9.5 

-2.2 

-5.4 

-24.5 


V-21 


SIGMA LEVEL 


COEFPICIENT/RUN 



PRESSURE-SMOOTHING COEFFICIENT 


:D = 0.0015 
(RUN SF-6) 


:D = 0.003 
(RUN SF-8) 


o 

• 

o 

-21.2 

1 

1 

1 

-14.4 

”27.9 



CM 

• 

1 

to 

• 

;.o 

i 

‘ 1 
1 

-43.5 







TABLE V-23: STORM INTENSITY VARIATIONS 


ITEM 

CENTRAL PRESSURE 

VARIATION (MBS) 


AVERAGE 

1 MAXIMUM 

Temporal Filtering 
(on-off) 

1 

0 

1 

0 i 

! 

Pressure Smoothing 
(2X variation) 

j 2.5 

4 

Surface Friction 
1 (2X variation) 

1 

1.5 

1 

4 ; 

Precipitation 

(on-off) 

I 3.0 

5 

Momentum Diffusion 
(lOX variation) 

3.0 

8 

Temperature Diffusion ; 

(lOX variation) 

1.0 

3 

i 

i Mountains (in-out) i 

1 i 

1.1 

5 

j Moisture Specification 
i (50^ variation) 

2.6 

6 

1 

1 Wind Specification 
(anal-balance) 

1-2 

2 ; 
1 

1 

; Radiation (on-off) 

j 

0 

0 

i Horizontal Resolution 

1 (3X variation) 

1 

6.3 

11 

] 

1 Vertical Resolution 
I (2X variation) 

0,5 

2 


VI . SUMMARY 


1 


A. Description of Modifications 

Sensitivity tests using model version PECHCV ( 63 x 63 
> 5 levels) have been completed. Not enough computer time 
was available to permit the series to be re— run using model 
version PECHPV (63 x 63 x 10 levels) . Modifications to the 
standard model, together with identifying run numbers, are 
indicated below; 

• Temperature diffusion coefficient (K^) 

= 10^ (standard) ; Run SF-3 
^ = 10^* Run SP-4 

• Momentum diffusion coefficient (Kj^) 

j — ~ (standard) ; Run SF-3 

— K,, = 10^; Run SF-2 

M 

I — ^M “ SF-1 

t 

I « Surface drag coefficient (C^j) 

r = 0.0015 (standard); Run SF-3 
D 

— Cp = 0.003; Run SF-7 
|i — Cp = 0.00015; Run SF-9 


I 


i 


VI-1 


• Pressure-smoothing coefficient (e) 

— e = 0.01 (standard); Run SF-3 

— e =0.02; Run SP-6 

— e =0.0; Run SP-5 

— e = f(t); Run SF-12 

• Precipitation (large-scale plus convection) 

— On (standard) ; Run SF-3 
— Off; Run SF-14 

• Temporal filtering 

— On (standard) ; Run SF-3 
— Off; Run SF-15 


• Multiple modifications 
Standard; Run SF-3 

No diffusion, pressure smoothing, friction, radiation, 
condensation, evaporation, or sensible heating; 

Run SP-11 

No radiation, evaporative/sensible exchcnges at 
surface; Run SF-20 


No diffusion (K^=Kj^j=0 . 0) ; Run SF-10 


Initial wind specification 


Analyr.ed (standard) ; Run SF-3 

Derived from balance equation; Run SF-16 

Initial moisture specification 

Parameterized using vorticity patterns (standard); 
Run SF-3 

— RH = 25%; Run SF-17 

— RH = 75%; Run SF-18 


‘ Tendency truncation/Terrain tapering 

Tendencies and terrain vanish south of 5°N, and 
increase SDf.-.j)thly with latitude to full value at 
20 “N (standard) ; Run SF-3 

Tendencies and terrain reduced by 50% south of 5°N, 

with latitudinal increases to full value by 20°N; 
Run SF-19 


Terrain specification 

Tapered south of 20°N to 5°N equatorward of which it 
vanishes (standard) ; Run SF-3 

Smooth earth; Run SF-13 
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B. Findings 


These sensitivity studies have demonstrated that the 
effects of each process, procedure and parameter contributing 
to the computer forecast need to be understood, ordered and 
exploited. Many of the choices are model depenaent. 

Physical effects are often small when compared to other 
sources of (possible) variation in a forecast. In fact, some 
computational/cosmetic devices can m^ (or overwhelm) certain 
physical effects. Consider, for example, pressure smoothing 
which has a greater impact on the forecast than surface friction, 
ountains or vertical resolution. Momentum diffusion impacts 
on the forecast as much as (modeled) precipitation. The 
coefficients for both pressure smoothing and momentum diffusion 

should be specified only after careful evaluation of numerous 
forecasts . 

Table VI-1 summarizes some of the variations between two 
forecasts due to specified modifications in a procedure, 
coefficient or effect. Precipitation amounts appear to be 
affected mainly by changes in the horizontal resolution and/or 
in the initial moisture specification. The average storm 
intensity is affected by horizontal resolution, too. Vertical 
resoj.ution increase (from five to ten layers) had minimal 


TABLE VI 1 : VARIATIONS IN FORECASTS FROM SPECIFIED MODEL CHANGES . 


r— — — 

ITEM 

DIFFERENCE BETWEEN TWO FORECASTS 

PRECIPITATION 
AMOUNT (%) 

AVERAGE STORM 
INTENSITY (MBS) 

KINETIC ENERGY | 
CHANGE AT LOV4EST 
LEVEL (%) 

RADIATION (ON-OFF) 

1 

0 

0.0 

TEMPORAL FILTER (ON-OFF) 

1-2 

0 

0.1 

FRICTION (2X CHANGE) 

2-3 

1.5 

14.6 

DIFFUSION (lOX CHANGE) 




- MOMENTUM 

5 

3.0 

66.9 

- TEMPERATURE 

3-4 

1.0 

4.7 

PRESSURE SMOOTHING (2X CHANGE) 

j 

4-5 

2.5 

14.6 

INITIAL WINDS (ANALYZED-DERIVED) 

5 

1.2 

12.7 

MOUNTAINS (IN-OUT) 

6 

1.1 

21.1 

PRECIPITATION (ON-OFF) 

N/A 

3.0 

14.0 

VERTICAL RESOLUTION (2X CHANGE) 

10-15 

0.5 

3.5 

HORIZONTAL RESOLUTION (3X CHANGE) 

80-100 

6.3 

57.7 

INITIAL MOISTURES (75%-25%) 

50-75 

2.6 

1 

19.2 



impact on storm centers. Pressure smoothing and momentum 
diffusion have as much effect on lows as precipitation or 
surface friction. With respect to the kinetic energy ' change 
predicted at the lowest model level, note that momentum 
diffusion has greater impact than tripling the horizontal 
resolution. Pressure smoothing has the same effect as does 
surface friction or precipitation. Both radiation and tem- 
poral filtering appeared to have only minor (if any) impact 


on a one-day forecast. 

The effects of the model modifications made were selectiv, 
with respect to model level, to features in the flow patterns, 
and to the type of condition being examined. The energetics 
can be affected, for example, without discernible differences 
being produced in the overall verification statistics. Many of 
the modifications produced effects which attenuated with alti- 
tude (precipitation mechanism; terrain; friction, pressure 
smoothing, horizontal resolution) . These tests were made using 
one data set. We would expect slightly different results If 
the sample was enlarged. Indeed, larger variations would bo 
e.xpected if the forecasts were of longer duration. Radiation, 
for example, would have greater impact on a one-week prediction 
some of the modifications were detrimental. This was expected. 

The purpose was to define a range of possibilities, rather than 
to improve the model's performance. 
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